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ABSTRACT 


The object of this thesis is to develop a procedure 
for calculating the dynamic displacement response of 
the platform of a fixed offshore structure acted upon by 
a regular wave train. The structure considered has four 
legs in a Square configuration, with the waves impinging 
normal to one side of the square. The procedure may be 
manipulated for use with other leg configurations and wave 
directions. 


The thesis is written in a manner useful to a designer 
of fixed offshore structures. An iterative procedure may 
be used to arrive at the critical wave-displacement 
combination. 


The types of waves considered are those for which Stokes’ 
3rd approximation applies. A modified version of Morison's 
theory is used for drag forces. Families of curves are pre- 
sented for use in predicting these forces and their centroids. 
The curves, plotted in dimensionless form, show the limits 
of applicability of Stokes! 3rd approximation. MacCamy- 
Fuchs diffraction theory is used for inertia forces. 


To calculate the dynamic response of the platform the 
authors fit the wibration problem of the structure to the 
classical theory for a linear, single-degree-of-freedom 
system. The wave forces are expressed in terms of a Fourier 
series. Displacements are calculated for twelve positions 
in one wave cycle, and a displacement curve is drawn. A 
simple tabular form is presented for calculating these dis- 
placements. 


Computer programs (in FORTRAN) are given for performing 
the bulk of the force and displacement calculations men- 
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tioned above. The theoretical procedure for unbraced and 
braced leg configurations is tested by experiment for 32 
cases. For these tests the ratio of maximum theoretical 
to experimental displacements varies from about 0.5 to 1.7. 


It is coneluded that the designer should incorporate 
a factor of ignorance of 3 in his design. 


The authors show that the largest waves are not 
necessarily the most critical for the structure. Rather, 
the critical design wave is a function of the frequency 
ratio (wave frequency/structure natural frequency). The 
designer must consider a range of significant waves which 
give frequency ratios of about 1.1 or less. 


It is recommended that calculations be performed for 
some real structures in ocean waves. From a comparison 
of the results with observed data a revised factor of ig- 
norance might be obtained. 


In Appendix I an approach to the solution for vortex 
shedding forces and frequencies in a train of surface 
waves is presented. Existing vortex theories are modified 
to apply to a vertical cylinder in finite amplitude waves, 
and boundaries are set on the predicted maximum forces. 
Experiments are conducted from which vortex forces are seen 
to fall within the boundary limits set by theory. The ex- 
periments also validate the proposed theory that the shedding 
frequency is a predictable integral multiple of wave fre- 
quency. It is recommended that a procedure be developed for 
determining platform displacements due to these forces. 


Thesis Supervisor: Donald R.F. Harleman 


Title: Associate Professor of Hydraulics 
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LIL. INTRODUCTION 


The problem of wave-induced vibrations in fixed 
offsnore structures suddenly was cast into the public 
news focus when on January 15, 1961, the United States 
Air Force Texas Tower Number Four collapsed 80 miles 
off the New Jersey coast in a heavy storm with the at- 
tendant loss of 28 lives[1). 

Prior to that incident, and since, no one has 
computed the dynamic displacements of such a structure 
acted upon by wave forces, to the best of the authors! 
knowledge. The design approach has been to make the 
natural frequency of the structure relatively high 
(about 50 cycles per minute) and then simply treat the 
wave force as if it were statically applied. This 
method is crude at best, and can be very misleading. 

During the winter of 1958-59 a platform motion 
study [2 | was made of Texas Tower No. 4. It was ob- 
served that the tower platform displacements were 
greater for 10 and 11 ft. waves (about +3 in.) than 
for 30 ft. waves. Obviously the highest wave was not 


the most critical for this structure. 





Y Numbers in brackets refer to literature citations, 
Appendix J. 
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The Air Force was interested in platform motions 
of Texas Tower No. 4 because they sacle to know what 
errors might be introduced into the installed aircraft- 
tracking radar system. (For this particular situation 
it was felt that rotational motion would probably be 
more of a problem than translational motion. ) 

This thesis develops a method of calculating the 
dynamic translation displacements of a fixed offshore 
structure in a regular wave train. The theory developed 
nas been experimentally checked by the authors. The 
experimental model was patterned somewhat similar to the 
United States Coast Guard's Buzzards Bay Light Station, 
located off the coast of Massachusetts eae 

The Coast Guard has initiated a long range pro- 
gram of replacing lightships by fixed offshore structures. 
The two completed to date (Buzzards Bay and Brenton Reef 
Light Stations) are braced structures with four legs. 
The U.S. Navy built a similar but considerably larger 
structure at Argus Island, near Bermuda. The U.S. Air 
Force built three Texas Towers, which were three-legged 
structures. Many other fixed offshore structures have 
been built in the past also, particularly in connection 
with drilling and mining operations in the Gulf of Mexico. 
The potential exists [4 | for many more offshore oil- 
drilling rigs to be built around the world in the years 


to come, 
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There has been at least one case reported [5] 
in which vortex shedding was important in waves. In 
that case after a one week storm a two foot diameter 
pile suffered a fatigue failure caused by transverse 
vibrations of 2.5 seconds period in 12 feet high 
waves of 13 seconds period. The authors of this thesis 
suspected that vortex shedding forces might have a 
noticeable effect on offshore structures. Some investi- 
gation of this subject has been included in this thesis. 
The authors concluded from their literature survey that 
no previous work has been done on the question of vortex 
shedding forces on a stationary vertical cylinder ina 


wave train. 








II. PROCEDURE 


A. Theoretical Procedure 


lig’ Hydrodynamics 


General 


This section deals only with those forces on a 
cylindrical pile acting in a direction parallel to that 
of wave propagation. These forces are termed longitudinal 
forces to distinguish them from the transverse or "lift" 
forces caused by vortex shedding. The transverse forces, 
although requiring consideration in a practical structure 
design, caused such small displacements of the model under 
study that it was not feasible to obtain experimental 
verification of an analytic approach. Because of this, 
the transverse force investigation, both theoretical and 
experimental, is included as Appendix I and is not dis- 
cussed further in this section. 

The ultimate design procedure would be to make use 
of a statistical analysis of ocean wave and energy 
Spectra in the development of the input forcing functions 


for the structure displacement problem. 
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Steps have been taken in this direction, for example 
[6], although the problems encountered, not only in ob- 
taining a sufficient scope of reliable data / 7 |; but also 
in applying this data to the development of a useful de- 
Sign criteria, have yet to be solved. For this reason 
the authors are restricting the investigation to a 
regular train of surface waves whose forces are capable 
of relatively simple mathematical Superposition, 

In the past a typical design procedure has been to 
decide upon a "design wave," which is the largest wave 
to be expected (within specified confidence limits) over 
the planned lifetime of the structure. Using applicable 
wave theory, one may then determine the wave forces and 
moments acting over an incremental length. These may be 
Summed up in a tabular form to give the total force and 
moment. [ 8, 3 ). Some designs, also, have considered 
a breaking wave of size somewhat smaller than the design 
wave [9, 10]. 

If the dynamics of the problem are considered it 
is obvious that this design wave may not be the wave 
which causes the greatest excursions (or stresses) in 
the structure. Rather it is to be expected that a wave 
(of significent height) whose fundamental or some low 
order harmonic corresponds closely to the natural fre- 
quency of the structure may be the critical wave from 


a design standpoint. This possibility has been indicated 
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where, over an extended period of observation, a 7.5 foot 
wave caused greater platform motion of Texas Tower No. 2 
than much larger waves did. [9] A similar observation 
is cited in the Introduction for Texas Tower No. 4. 

A serious situation also may exist when the wave 
length is approximately equal to the leg spacing, for 
in this case the fundamental wave forces (and possibly 
some harmonics) will be in phase. 

The proposed procedure, therefore, presents a 
method whereby several waves of various lengths and 
steepness may be considered over the range of possible 
water depths. The critical displacement-wave combination 
is then arrived at by an iterative process. To allow for 
a solution by this method within a reasonable length of 
time families of curves are developed such that total 
force and moment may be determined by a simple muiti- 
plication. This is intended to supplant the tedious in- 
cremental tabulation and summation mentioned previously. 
The curves (Figures III, IV, V) are presented in terms 
of dimensionless parameters such that they may be applied 
to any situation for which the proposed theory is appli- 
cable. 

A pile spacing greater than 10 pile diameters has 
been found exp@rimentally to be large enough sucn that 
proximity effects are negligible. /11/] This spacing 


is exceeded, not only in the vast majority of offshore 





structures of the type considered, but also for the 
model which was the object of the authors! experi- 


mental study. 


Wave Theory 


H = wave heignt, measured from 
crest to trough. 


L = wave length. 
d = water depth. 


We height of wave crest above 
stiil water level. 


Several theories on wave motions have been developed 
and tested experimentally with varying degrees of 
success. These theories are found to be applicable 
over finite ranges of wave parameters and water depths, 
with no one theory satisfying all conditions. Since 
the wavés of interest in this study are of the re- 
latively shallow water finite amplitude type, two 
theories seem to be of particular interest - namely, 
Stokes! third approximation and the Solitary wave theory. 
These theories are discussed in detail elsewhere (for 
example [4, 12, 13, 14)) and will not be repeated here. 
However, for purposes of evaluating this study, the 
following limits of applicability may be noted from 


these references. 





Stokes' 3rd approximation: 


a.) (9 fH) 4 0.625 
b.) (H/L) 0.142 


IN 


Solitary wave: 


a.) waves at or near breaking point 


b.) (d/l) —~> 0 
Tneoretical limits on ocean wave parameters: 


(7 fH) & 0.75 
(1/7) A= no), WANe 


bee) 
<i 

ec.) H = 0.78 (a-(H- 7% )) for shallow water 
sia) 


b 


max 


d.) from the above, (H/d) 4 0.653 


Using the maximum value of (7 fH) from Stokes! 
theory and the value of ax above, one sees that 
rd) is 0.603 by Stokes' theory. There is a 
gap between this value and 0.653 for which the soli- 
tary wave theory would appear to be more applicable. 
However, waves have been observed breaking on Martha's 
Vineyard with values of H/d = 0.588. [15] 

Because this present investigation must be limited, 
it was decided to use only Stokes! 3rd approximation. For 
the solitary wave and other theories the method of approach 


in determining forces and force centroids basically would 





be the same, with the curve families of this thesis en- 
larged to include their particular ranges. Bretschneider 
[14) has developed graphs which appear to be a good start- 
ing point for determining the range of applicability of 
the several theories, as well as values of H, I 
maximum crest velocity (Una? for design waves. These 
values could serve as first trial inputs for the pro- 
cedure to be presented here. For a conservative estimate 
of the force exerted by a breaking wave oan a piling a 
value of 4 times the force as determined from Stokes' 
theory for a wave of the same height may be used. [16 } 
For reference purposes the applicable formulas 
from Stokes! 4rd approximation are included here in 
dimensionless form. [17] The reader is referred to 
the list of symbols at the beginning of this thesis and 


Figure I for use with these equations. 


Let C = wave celerity, ft./sec. 
J = 2nd/L (2) 
The vertical displacement of a particle from 
its mean position is given by (2) 


5 ca [SRE YD] ono [8 ty SR SED on 
(2) 
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Figure I. Definition Sketch of a Wave 


At the wave crest (2) becomes 


Ymax i sinh(J S/d) _ 3 ww sinn(2J S/a 
a. “> aH a “sinh J * 8 oat ah aoe (3) 


(sinh J) 
and at the trough, 
Vtrou ne ut -Sinh (J S/d) 4 3 sinh (2J age (4) 
ie = 2 d/H sinh J 8 L/H (euine. 


At the still water level (S/d = 1), (2) gives the wave 


profile 
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Lm raya | sin 9 -|3 Vii SE Pe, cos bog (5) 


(sinh J) 


From (5), with 9 = 90° the height of the crest above the 
still water level is 


8 
(sinh J) 


a stall +h oa mea af? (6) 


and for 9 = 270° the distance to the trough is 


=H 
Lo - Saf +8 op inh (23) ; (7) 





rs) 
(sinh J) 


The wave celerity is 








8(sinn J) 


The particle velocity at any depth S/d is 


u_ or j eeu S/d) /esn 9 -| —* 1 (al 
yaa WH Vea sinh J I/H"(sinn 3)* 


)] cos 20 + a3] ° 29 - if 


sinh J) 


At the wave crest (9) becomes 


at 


(9) 





2 2 
oo. 2 i Ls 2(cosh 2J)"+ 2cosh 23+ 
Vv gd Vd tenn 3 $2 + gL i (8) 


cosh (2J S/d 


(sinn J)* 


rr om |e, 





wMAX ul ~ a <a od 1 (- _ cosh (2J ) 
/ ga WE aa bailed WF (sinh J) 2 (sinh J)” ~ 


ae osh(2J S/d } 
: rae | ae 3) I a : 


At the wave trough (9) becomes 


Ww 
trough mw ¢  ( cosh(J S/a) , 1 _ 1,3 cosh(23 S/4)) 
L/H’ } sinh J L/H° 2\" 2 2 | 
inh J) (sinh J) 
Y, gd Jed ie | 
nw 1 7 cosh(2J S/d & 
L/H*2 (sinh J) J 


Force Theory 

Using Stokes! 3rd approximation for wave motions, 
D.R.F. Harleman and W.C. Shapiro of the Massachusetts In- 
stitute of Technology Hydrodynamics Laboratory have de- 
veloped a force theory which shows good agreement with ex- 
periment. [ 17) In developing this theory they applied 
certain modifications to Morison's Force Theory for drag 
and the MacCamy-Fuchs Diffraction Theory for inertia forces, 
In this thesis the authors will utilize the approach as pre- 
sented by Harleman and Shapiro. Families of curves will be 
developed from their equations, these curves to be used for 


the iterative design procedure mentioned previously. 
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As will be shown the only forces of interest are 


the maximum values of the drag and inertia components for 


both crest and trough regions. 


where 


maximum drag force in crest region 
maximum drag force Intrough region 
maximum inertia force in crest region 
maximum inertia force intrough region 


cylinder diameter 


2C . 
~7{D a° “D trough B 
(Fp), 2 


14/4 42 
= J a ae a(S'/d) 


1b 


d 2 
Uu 
B= f een A a(s'/a) 
0 


(12) 


(13) 


(14) 


(15) 


and S' is the vertical distance to the actual particle 


position and is therefore equal to the mean particle 


position S plus or minus the vertical displacement y. 


1) 





ee re 


-_ = —— MD az 





For the crest S'amtS tye. 2 (16) 


and for the trough S'«S - gous (a7 ) 


It can be seen that (14) and (15) do not lend themselves 
to regular integration, but must be solved graphically or 
by some other means. The steady state drag coefficients 


and C to be used in (12) and (13), respectively, 


Cp D trough 
are determined as functions of Reynolds numbers, defined 


as follows: 


2 
a LS». BYR 08 


} ut 
IR y_ trough p ot £2 2 (19) 


trough — ) . (Lo 


nee. [3 
where Cac and u trough are the root mean squares 


of the velocity distributions for the crest and trough, 


respectively. 


|—~_ P) 
2 Aga 
WF asi a ; = 


2 
5 Bed 
J Utroueh a+ (y-H) -_ 


In (18) and (19), A and B are determined from (14) 
-H 
and (15), respectively, and 7fo» = from (6) and (7), 
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Figure II. Steady state drag coefficient versus 
Reynolds number for circular cylinders, 


(eee 10) 


respectively. Figure II is a plot of steady state C, vs. IR 


for circular cylinders. 
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F, = LB H sinn I cos ¢ (22) 
T p* H ; 
(F.), = Ot sinn I cos @ (23) 


wnere I x= J + DA sin 9 (24) 
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Figure II. Steady state drag coefficient versus 
Reynolds number for circular cylinders, 


pln 10). 


respectively. Figure II is a plot of steady state C, vs. IR 


for circular cylinders. 


2 

— WB a H sinn I cos ¢ (22) 
T De H : 

(F,), = at ta sinn I cos @ (23) 


where I = J + DA sin 9 (24) 


F 
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and 9 is that angle of wave position (0° <= 9< 180° 


3 6 
for (22), and 180 <0 < 360 for (23)) such that F. and 


ie 
(Fr), are maximized, respectively. Because 9 will be 
very nearly equal to 0° and 180° for these situations, 
it can be seen from (24) that little error (and great 
simplification) is introduced by letting 9 = 0° and 

180° in (22) and (Cor; respectively. 

T p° H 


This gives FL = -(Fy), = tanh J (25) 


Let S' = center of action of the force F, and 





D 
(S'), = center of action of the force (FL )t 
Nr d 
S! "ef Umax” 
Then = J S'/d “gd a(S'/d) (26 ) 
“| J 
A 
-H 
(S') TT 2 
ee trough 
=ar ea /a) (27) 
° 
v5) 
8, = center of action of force F, 


2 
al) Wek a | 
5, > oF cosh J =o s/o 1+ J sinh J cosh J (28) 
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Harleman and Shapiro have shown that S' can be 
assumed constant for O +6 © 180°, and (S'), can be 
assumed constant for 180° = 9 = 360° with little loss in 
accuracy. Also, with the assumption that Sy is equal 
for crest and trough by (28), it is reasonable to assume 
that 1t is constant for any value of @. 

By use of (5,7,10,11,14,15,26,27) values of A and 
B for (12,13,18,19) anda _ and a for (26) and (27) 
were computed for families of the parameters H/L and d/L 
on an IBM 7090 digital computer. 

The results are plotted in Figures III, IV, and V 
in dimensionless form. Values of H/d equal to 0,603, 
which was mentioned as the approximate limit of applica- 
bility of Stokes' theory, and 0.653, which is the maximum 
value this parameter may obtain (page 8), are also 
Shown in these figures. These values indicate the boundary 
for use of the figures, as well as the boundary up to which 
. the solitary wave theory could be used to extend the range 
of applicability. The computer prerram, representative 
output, and comments are included as Appendix D. In ad- 
Gitien te A, B, ., and Ca the following are given 


as outputs: 
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H/a = 9,653 








7am 0.6035 (limit 
of Stokes! theory ) 


H/d = 0,653 (theoretical 
limit ) 


(s" eS (4) | vie Itt 
0.2 Drag Force Centroids for Wave Height/Wave Length 
| (H/L) and Water Depth/Wave Length (d/L) Ratios. 
(8 /a) « centroid/depth ratio for crest | 
3 'ercurs /a) = centroid/depth ratio for trough 


O.1 


| | a/L | don 4/13/62 | 
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Velocity profiles are also given for mex , *max 
_ =" 


“trough . A Yt rough® i, SO 
yea ~ ae 


Design Force Procedure 
It has been shown [aaa that drag force is a 
function of sin aon while, for the assumption that Fr 
occurs at 9 = 0°, 180°, inertia force is a function 
ef cos 9. The forcing functions for use in the dif- 
ferential equations of platform displacements will there- 
fore be of the following form: 


aa 6 = 180° > 


E. F = Fp sin “9 + F,; cos 0 (29) 


180° = @ * 360° 


F = (Fs sin“@ + F. cos O° (30) 


6 


For convenience the formulas and procedures used 
in the determination of these forces-and their points of 


application are summarized below. 


on 





<waiZ Bam “pA | 
“7 2 (12) 
(F.), = _y D a? SD tro h B (13) 
Dt 2 
T p* is 
F. = it! tanh 3] (25) 
Cy from Figure II for R = . S60 3 . (cy 


a‘ + 7/2 


D gdbB 


CD tr ough from Figure II for IR = = %> 
— 
A from Figure IV 
B from Figure V 
Mo 1 j 3 mw sinh (2d 
= L + (6 ) 
ad 2 d/H 8 L/H (sinh J) ; 


or from computer output ade plotted) 





-H 
( sao - 5 (31) 
S! (S!# )t 
a7 and q from Figure III 


Sais See 2) (28) 


oe 


Ze 


(19) 





2. Vibrations 


General 

Once the wave forces on the legs have been deter- 
mined the next problem is to calculate the resulting 
displacement of the platform as a function of time, or 
wave position 9 with respect to the front (upstream) 
legs of the structure. 

This thesis considers only the longitudinal dis- 
placement x of the platform, taken as positive in the 
direction of wave propagation. It has been observed 
that the transverse vibrations in the experimental runs 
were relatively negligible, with an amplitude of perhaps 
three percent of that of the longitudinal vibrations: 
For the dimensions of real structures the transverse 
vibrations may be as important as the longitudinal 
vibrations. 

This thesis does not consider wind force and its 
effect on platform displacement. Certainly with large 
waves there will also be high wind, and the wind-caused 
displacement must be added to that induced by the waves. 
In a steady wind the problem is one of statics, not 
vibrations, unless vortex shedding is important. A useful 
discussion of wind loads is given in [4] . Water currents 


can be treated in the same way. 
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The structure considered has a square platform of 
weight W, supported on four (N=#4) legs of length X and 
diameter D, arranged in a square configuration of leg 
Spacing h on a side. The case of no leg bracing is con- 
sidered first, after which the case of added bracing is 
considered. 

The investigation is limited to the case of the 
waves impinging normal to one face of a Square con- 
figuration of legs. Once the reader understands the 
concept of this procedure he should be able to manipulate 
the procedure to fit any practical leg configuration and 
a limited number of angles of wave impingement. When the 
direction of wave propagation is not parallel to an axis 
of symmetry of the braced leg configuration, the calcu- 
lations become complicated by the facts that generally 
the platform translation is not in the direction of the 
applied force and that rotation of the platform may be- 
come Significant. 

The legs of the thesis model are fixed at the base. 
This condition was assumed for the Argus Island structure 
| 197] . In the case of Texas Tower No. 4 the footings 
beneath each leg were 25 feet in diameter filled with 
cement, and sunk or embedded into the ocean floor to a 
depth of 18 feet. After the collapse divers reported 
that the footings were apparently in 06a condition with- 
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out any evidence of fracture, movement, or scour [20] : 

The opinion was that the footings were not cocked out of 
a perpendicular plane with the ocean floor [ 21] P In 

other words, the assumption of fixed footings was valid 

for Texas Tower No. 4. 

The legs of the thesis model are fixed at the plat- 
form. This assumption was made in the design of Texas 
Tower No. 4. Later investigation [2] indicated that the 
actual restraint condition was closer to being pin-ended 
than fixed. 

As a practical matter in the consideration of 
vibrations, it should be pointed out that large quantities 
of surging oil or water in tanks may amplify tower motion. 
This problem generally can be reduced by installing 
several small tanks in lieu of one large tank and by in- 
Stalling internal baffling in the tanks. 

Reference { 20 ] on the collapse of Texas Tower No. 4 
provides much useful and interesting food-for-thought 
for the designer of fixed offshore structures. 

The reader interested in experiments with models 
Should realize that it is not possible to scale model 
results to fit the prototype. The problem is that drag 
and inertia components of the wave force do not scale 
in the same manner. Consequently the reader should 


center his interest in the concept of this thesis, and 
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not on the numerical results for any particular 


run. 


Classical Vibration Theory 


The approach considered will be to fit the 
vibration problem of the structure to the classical 
theory for a linear, single-degree-of-freedom system. 
There are many sources in the literature where one 
may find this theory. The authors have chosen to 
follow the notation of Den Hartog | 22]. 

From Newton's Law one can arrive at the follow- 


ing differential equation of forces: 


mx + Ox + kx = Po sin Wt (32) 


where displacement, a function of time 


gn/ ot = velocity 

x/at© = acceleration 

mass with single degree of freedom 
damping coefficient 

spring constant 

magnitude of applied force 

circular frequency of applied force 
time 


ct fo” ms a ie 


The steady state solution to (32) can be written as: 


x= x, sin(wt - Z) EB) 


where x, and % are as yet unlmown, and: 
{ 
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Xo = displacement amplitude 
@ = phase angle by which x lags Py sin @Mt 


Insert (33) into (32) to obtain: 


Oo 
2 2,° 
(CW) ” + (k - mw”) 
@ = arc tan ee (35) 
-~ m 


From consideration of free vibration without damp- 


ing one can find that the undamped natural frequency 


Die |= (36) 


From consideration of free vibration with damping 


W.. is: 


‘aT 


one can define critical damping Cc. as: 


C,=2/ mk = 2mw, (37) 


For values of damping C less than C. the free 
vibration is truly oscillatory. For values of C 
greater than C, the free "vibration" is not oscillatory, 
but rather the mass creeps back exponentially to its 


equilibrium position from an initial displacement. 


as 


ay 





Insert (36) and (37) into (34) and (35) to obtain: 


Mao. (38) 


LL —- — 
XSTATIC 2 d 2 
h- a + | 
n Cc 


Po 
where Xonartc * iq 


» 


(39) 


(40) 


It will prove useful to have numerical tables for 
expressions (38) and (40) for entering arguments of 2 
and small values of C/C.- Table I is presented for 

x 


values of —~———-._ Table II is presented for values of @. 
*STATIC 


Unbraced Structure 


Mass m must represent the mass of the platform 
and a certain portion of the mass of the legs. By 
applying Rayleigh's Energy Method to this case (Appendix 


A) it is possible to compute: 


“— = (Ww + sot) (421) 
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TABLE I - MAGNIFICATION FACTORS (38) 
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TABLE II - PHASE 
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MH Oo 


OCOOVDODOKFNN FAONDNLUNDRFKRH EF IO 
WN 
WW 


2.12 


OM Ut fdas 


OOO0O0OF 
MUI DMO O 





where ¢ 


gravitational constant, normally 
32.2 ft/sec.2 

weight/unit length of leg 

number of legs 


az 


If the configuration had been such that the legs 
were considered to be pinned to the platform the 
fraction to be used in the second term of (41) would 
be 33/140 (Appendix A). 

It is shown in Appendix A that the added mass 
effect (of accelerated water mass) is negligible in 
still water. For forced vibrations there is the further 
consideration that the water particles are moving with 
wave orbital velocities. For this latter situation it 
is simply assumed that added mass is negligible. Later 
in the experimental results (Appendix G) it will be 
seen that this assumption is valid. 

For damping C it is necessary to resort to experi- 
ment to determine values of = for the plastic model of 
this thesis. The damping Someidie of structural damping 
in the plastic and viscous damping by the water. Again 
from Den Hartog | 22 |for the case of free vibrations with 
viscous damping one finds that the logarithmic decrement 
JS fits the following relationships: 


a 
a Log, ro (42) 
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where x, and x, are two successive positive (or 


negative) decaying displacement amplitudes. 


a 


C 
ie (43) 
Ce ‘ 





For values of c/c. up to 8%, and somewhat larger, 
the following approximation is "exact" within sliderule 


accuracy: 


fe aed (4424) 

C 
With this expression one can calculate values of 
structural damping by observing free vibrations in air, 

assuming that the air causes negligible damping. 


For completeness at this time it is noted that: 


(45) 





where q = damped natural frequency 


With this expression one can calculate experimental 
values of Wo: 
It is shown in Appendix A that in still water for 


these particular experiments the theoretical viscous 
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damping is 11% maximum of the structural damping. In 
the steady state vibrations there is the further con- 
sideration that the velocity associated with viscous 
damping is not x but is the velocity relative to the 
water particles, which are moving with wave orbital 
velocities. For these reasons viscous damping by the 
water is neglected. The experimental results (Appendix 
F) support this action. 

In Appendix A the following expression is derived 


for spring constant k: 


k = s (46) 
where E = Young's modulus 
and I = moment of inertia. associated 


with bending of one leg 


Actually k must be determined experimentally for 
the plastic model of this thesis because a unique 
value of E cannot be obtained from a handbook for the 
plastic. 

From (36,41,46): 


i. 12gNEL (47) 
1 
n pr + se. 4) 


De 





The forcing function in (32) is oa sin wt. 

The expression for the wave force is not this Simple, 
but if it can be represented by a series of terms 

of the form Po 8inwt one can find an equal number 
of terms for x, each with its own phase angle. By 
superposition of these terms one can calculate the 
platform displacement x. 

The forces acting on each leg are equal in magni- 
tude. The forces on the front legs are in phase 
(because of the direction of wave impingement on the 
structure). The forces on the rear legs are likewise 
in phage with each other, but they lag the forces on 
the front legs by a lag angle f,,, which one can deduce 
to be: 





gat ~. 220 h degrees ~ ca radians (48 ) 


For the drag force acting on one leg it has been 


stated (29,30) that in the range 0° 2 9 & 180° 


2 


the force can be expressed as F, sin°-Wt, At 9 = 90° 


D 
the line of action of this force is at a height s! above 


the bottom. For 180°42 9 < 360° the arag force is 


2 = 
(For ough Sin” Ww t, which acts at a height (3'),- The 


inertia force acting on one leg can be expressed as Fy 


cos Wt, acting at a height 8, when 9 = 0°, 180°. 


Be, 
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Only a fraction of these forces can be considered 
to act directly on the platform. Let PG represent 
the force required at the level of the platform to 
produce the same static deflection as force P acting 
at height a on a leg of length 2. An influence 
fraction * o/ has been computed in Appendix A, with 
the numerical results given in Figure VI. 
It has been stated earlier that F, can be con- 
sidered to act at s' throughout the range 0° = 9 € 180°. 


This means that in this range: 


3! 
| h 
y (49) 


> | 


Use Figure VI, and write 


P 
F) = Fy-p~ (50) 


With similar reasoning: 
a 


ee 18 a positive number. ) 
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FIGURE VI - INFLUENCE FRACTION 
mF a ele Sap omOr 


‘Praction 52 of Force P to be.Considered Acting On 
Platform When P Acts at Height.a on an Unbraced 
Leg of Length Q. 
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pS (53) 
F, = F x (54) 
ale SS 


Now one can write expressions for the force acting 
on the platform due to the force, acting on one front 


leg. 


1 
force = “Fin t + F, cos w t (180° wt 4 360°) (56) 


force = F. sin-wt + F,coswWt (O%wt 4 180°) (55) 


_& Fourler series is calculated in Appendix A to re- 
present this force. For the cases considered in this 
thesis only five terms are significant. The following 


expression results: 


F, - F 
F(t) = + + sa (F, + F,) sin Wt 
4 
: Taq (Fy + F,) sins Wt + Fz cos wt 
end 


2) cos2 Wt (57) 


The angle Bat, has been defined (48) for a harmonic 


force varying as Wt. One can deduce that if the force 
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varies as nWt the appropriate lag angle is 
nPrr: 


lag angle = nf, (58) 


At this point all of the necessary theory is avail- 
able for computing displacement x as a function of 6 
at the front legs. One has simply to set up a con- 
venient method for a hand numerical calculation. 
Fortunately a simple tabular method can be found. For 
purposes of illustration it seems best to show this 
method in conjunction with a sample calculation. The 
reader is therefore referred to Appendix B. 

The numerical calculation can also be performed 
on a digital computer. Such a program was used in the 
preparation of the Results, Section III. The program 
is given in Appendix D. 


Braced Structure 


The solution to the problem for a braced 
structure may bé approached by applying rough modifi- 
cations to the solution for the unbraced structure. 

For mass m it was found previously (41) by 
Rayleigh's Energy Method that 13/35 of the mass of the 


legs should be added to the mass of the platform. 


at 





Assume that the same ratio holds for the bracing and 


legs of the braced structure, Define the ratio ra as: 


_ mass of legs + mass of bracing 
im mass of legs (59) 


The assumption is: 


m = =(W + sé 0,N n wi ) (60) 





legNEL 






n 3 3 r 
N-W + sé Nw) 
In the design of the Argus Island structure the 
fraction E was used in lieu of 32 .[ 23 | It appears 


that this fraction was assumed, a derived. 
Damping for the plastic model must be determined 
experimentally, as before (44). 


"In most metal structures, nearly 

all the damping observed can be 
“attributed to the technique of 
fabrication. For example, a good 
weld can be distinguished from a 
poor one by its effect on the part's 
ringing after impact: the poorer the 
joint, tne more damping. Riveted 
joints generally contribute more 
damping than welded ones, and bolted 
joints more than riveted ones." * 


Hamme, R.N., [ 24, p. 37-2] 
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"Bare steel structures, particularly 
welded rather than riveted, have very 
low damping ranging from perhaps 1.0 
to 3.0 per cent of critical damping.” * 


If a designer is in doubt as to what damping exists 
in his structure he will generally be well advised to 
assume the smallest reasonable value in order to be on 
the “safe side." 

Stiffness k remains as expressed in (46) for the thesis 
model because no bracing was added in planes where it would 
have stiffened the structure for resistance to longitudinal 
forces. If such bracing had been added, experimental dis- 
placements would not have been large enough to be ob- 
servable. 

For a normal braced structure one must compute ap- 


propriate values for stiffness k and influence fraction 


P 


7 A good textbook on this matter is | 25 } » Probably 


two values of Po will be sufficient to fair a curve, be- 
cause the hy will be very similar in shape to that of 
Figure VI. 

The case of Texas Tower No. 4 provides an interesting 
lésson in engineering with regard to stiffness furnished 
by the bracing. On this tower the bracing consisted of 
pin connected tension members. It was eoncluded [ 2} 


that the underwater bracing was essentially ineffective 


. Housner, G.W.,{ 24, Pp. 50-20] 
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for platform excursions up to +3 inches because of the 
large clearances in the pin connections. The lesson is 
that joints should be welded, which is the method usually 
used. 

It is seen in (12,13) that drag force varies directly 
with leg diameter D. It is seen in (25) that inertia 
force varies directly as Dey For each of these forces 
the line of action is independent of D. Assume that if 
bracing existed whose long dimension was parallel to 
the direction of wave propagation the wave forces on 
this bracing would be negligible compared to the forces 
on the bracing whose long dimension was parallel to the 
Wave crests. Let X,, be the cumulative length of all of 
the hydrodynamically significant underwater bracing 
(below the still water level), of diameter D,. For 
Bimplicity assume that the wave forces on a length 
d(d = water depth) of bracing are the same as if the 
bracing were uniformly distributed in one additional leg 


of diameter D These assumptions lead to the following 


B° 
relationships for forces acting on the braced structure: 


(Fon = "p-Fp? (Pilg ™ "p-Fy> (alg ™ tp-Fo (62) 
where ry 3 ee 220 —) (63) 
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Melee Fy Folge tgede (64) 


D 
where —— oats Na >) (65) 
With these modified values for the forces to be 

used in (57) one can make a hand calculation for dis- 
placement x as in Appendix B, or a computer calculation 


as in Appendix D. 


Vibration Stresses 


Up to this point no mention has been made of 
stresses. It should be obvious that if platform dis- 
placement is greater for the dynamic situation than for 
the static situation the ettendant vibration stresses 
throughout the structure will Pino be greater by the 

“same ratio. (Deflections are "small," and the stress- 
strain relationship is linear elastic. ) 

The maximum dynamic displacement has been ob- 
tained from the plots of x vs @ that resulted from 
carrying out the foregoing procedure. AS a result of 
the calculation for stiffness k there is available a 
moment (and force) distribution (which is simply re- 


lated to a atress distribution by simple beam theory). 


The stress distribution is for a selected force, which 
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results in a calculated platform displacement. This 
static displacement is the denominator of the sought - 
for ratio. All values in the static stress aie oe - 
bution are multiplied by this ratio to arrive at the 
vibration stress distribution. 

In addition to the vibration stresses one must 
also take into account static stresses due to sueh 
eauses as wind forces, current forces, column loads, 
thermal effects, fabrication, ete. These stresses 
can be superposed to arrive at the total stress at any 


point. 


B. Experimental Procedure 


The experimental work of this thesis was carried 
out in the Wave Tank of the Hydrodynamics Laboratory 
of Massachusetts Institute of Teehnology in March 1962. 
The primary objective of the experiments was to 
measure platform displacement of a model of a fixed off- 
shore structure under the influence of a regular train 
of waves. Thirty-two runs were made. In addition, runs 
were made to measure the longitudinal and transverse 
wave forees on a single vertical cylinder. More will be 
gaid later about the various runs and the data that was 


collected. 
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Most of the equipment used, with the exception of 
the models, was in existence prior to these particular 
experiments and has been described in detail in 


several references |17,26,27,28 | . 


Wave Tank and Wave Generator 


Figure VII is a photograph of the wave tank, 
which is 90 feet long, 24 feet wide, and 3 feet deep. 
Basically it is a steel structure, but the walls and 
40 feet of the bottom are plate glass. 

At the near end of the tank there is a hinged 
plate wave —— hinged to the tank bottom. A 
rod connects the Gende end of the plate to an eceentric 
point on a rotating wheel. Wave height is governed by 
adjusting the eccentricity. Wave period is governed 
by adjusting the speed of rotation of the wheel. The 
speed can be determined by a photocell arrangement 
wired to an electronic counter. 

At the far end of the tank there is a 35 foot 
long beach arrangement, whose purpose is to absorb energy 
and prevent (or reduce) undesirable ware reflections. 

For these experiments, however, it was decided to avoid 
wave reflection altogether by starting the wave generator 


for each run and recording data before the refleeted wave 
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Figure VII - Photograph of Wave Tank 





Figure VIII - Photograph of Braced Structure Model 
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had time to return to the model. After some trial and . 
error with the waves to be used it was decided that data 
would be reeorded between the passing of the eighth 


and ninth wave crests at the test stand for all runs. 


Test Stand 


The test stand is seen in Figure VII. It is 
located approximately 40 feet downstream from the wave 
generator. It consists of a triangular frame of steel 
angles supported by three legs of 3-inch pipe straddling 
the tank. The stand is guyed to the floor by steel 
cables to inerease its rigidity. A steel angle was 
added to connect the stand to the laboratory wall to 
improve transverse rigidity, with a consequent improve- 


ment in the rigidity previously supplied by the cables. 


‘Models 


Figure VIII is a photograph of the-model of a 
braced structure. It was fabrieated of plexiglas by 
Forest Products, Inc., of Cambridge, Mass. 

The legs are 3 inch (actually 0.505 inch average) 
rod 42 inches long, with 16 inches on a side between 


leg centers. The leg diameter was decided upon as a 
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compromise between the requirments for a large diameter 
for adequate experimental forces and for a small diameter 
for adequate flexibility in order to observe deflection. 
The bracing is 4 inch rod, with horizontal members 

located 14 and 28 inches above the base. The platform 

is 5/8 inch thick (actually 0.610 inch average) by 
18-1/16 inches square. A 3/8 inch bolt in the center 

of the platform serves as a pointer for displacement 
measurements and holds the weights (to be discussed later). 
The base of the model is } inch thick by 30-1/8 inches 
by 22 inches, with beveled leading and trailing edges. 


The model also was used in the unbraced condition. 


Single Cylinders 


In order to measure longitudinal and transverse 
wave forces two different aluminum cylinders were used, 
They were attached to the portal gage by a flange ar- 
rangement, as seen in the photograph of Figure IX. The 
% inch cylinder cleared the tank bottom by 1/8 inch. 
The 1 inch cylinder cleared the tank bottom by 12-5/8 


inches. 
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Figure IX - Photograph of Single Cylinder in Position 


for Measuring Transverse Forces 
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Movie Equipment 


A 16 mm Bell and Howell movie camera was used 


to record platform motion and time, as seen in the 
photographs of Figures X and XI. Camera speeds of 


24 and 48 frames/second were used. 


Recorder 


A four-channel direct writing model 150 San- 
born recorder was used to record wave profiles and 
force traces. A paper speed of 25 mm/sec was found 
to be most satisfactory. Time marks could be made 
in the paper margin upon signal from a manual push- 


button. 


Wave Profile Measurement 


Wave profiles were measured by a resistance 
wave gage, which can be seen in Figures IX and X. 
The gage consists of two platinum wires (0.008 inch 
diameter) insulated from each other and oriented 


vertically with a spacing between them of 2 inch. 


When the wires are partially immersed in water an al- 


ternating current flows which is essentially directly 


proportional to the depth of immersion. The current 
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Figure X - Photograph of Braced Structure Model 


in Position for Measuring Displacement 





Figure XI - Photograph Showing Scheme for Recording 


Platform Displacement 
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produces a signal which is amplified and recorded by 
the Sanborn recorder. 

The gage can be raised and lowered by a graduated 
rack and pinion arrangement. It is calibrated statically 
by raising and lowering it through known distances in 
still water. 

Two gages are used to determine wave length. The 
gages are raised until their probes are immersed about 
& inch in the wave crest. The distance between the 
gages is adjusted by trial and error to one wave length, 
at which time two crest marks will be made simultaneously 


on the Sanborn recorder. 


Force Measurement 


Horizontal forces were measured by a portal 
gage, seen in Figure IX directly above the cylinder 
flange and oriented in this case to measure transverse 
forces. The gage measures horizontal force in terms of 
shear deflection. Il1ts design is such as to make the out- 
put essentially independent of the distance to the point 
of application of the force, or bending moment. A 
Schaevitz linear variable differential transformer is 
used to convert the shear deflection of the bottom plate 
of the gage into an electrical Signal, which is amplified 


and recorded by the Sanborn recorder. 
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The force gage can be calibrated in the positive 
and negative directions by the use of weights suspended 
on a string which is run over a pulley. Two pulleys 
are seen in position on the sides of the tank in 


Figure IX. 


Model Runs 


Water depth was maintained at 27 inches through- 
out the experiments. 
Four different waves were used, with the charac- 


teristics as listed in Table IIL. 


TABLE III 


Experimental Wave Characteristics 


Wave H (feet ) L(feet ) a7 H/L 


A 0.240 Daa 0.148 0.0158 
B Orel 10.86 CO IC g 0.0289 
C orool 8.82 OV295 0.0409 
D 0.404 {PS oe C07 ©0552 


The model was used in two configurations, un- 
braced and braced. For each configuration and wave com- 


bination four structural natural frequencies were in- 
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vestigated by adding weights to the platform. The 
weights used were 0,3,6, and 10 pounds. 

A total of 32 model runs was made. For each run 
platform displacement and time were recorded by the 
movie camera, and wave profile was recorded by the 
Sanborn recorder. A time mark was made on the re- 
corder paper at a known time as read on the clock that 


appeared in the movie. 


Single Cylinder Runs 


Single cylinder runs were made for each wave 
on a 4 inch cylinder to measure longitudinal and trans- 
verse forces. Runs were made for each wave on a 1 inch 
cylinder to obtain additional data on transverse forces. 
For each run wave profile and force trace were 


recorded. 


Miscellaneous Measurements 


Movies were taken of the free oscillations of 
the unbraced and braced structures in air and in water 
for all combinations of the four added weights. Dis- 


placement and time were recorded for three cycles. 
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An experimental value of spring constant k was 
determined. 


Experimental values of the influence fraction 


P 

p- (Figure VI) were determined as a function of 5. 
The weight of the plastic in the model was re- 

corded. 


Water temperature was recorded. 


20 





III. RESULTS 


Figure XII presents the theoretical and experi- 
mental results for platform longitudinal displacement x. 

The terminology "Run AO" indicates a run with 
wave A and O lb. of added weight on the platform of the 
unbraced model. The letter "Z" specifies the braced 
model. 

Experimental wave characteristics are given in 


Table IIitI. 


Values of structural damping _ used in the cai- 
c 


culations'are given in Table X. 
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Figure XIIb - Platform Displacements for Wave A 
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Figure XIIc - Platform Displacements for Wave A 
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Figure XIId - Platform Displacements for Wave A 
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Figure XIlIe - Platform Displacements for Wave B 
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Figure XIIg - Platform Displacements for Wave B 
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Figure XIIhn - Platform Displacements for Wave B 


OrnZ2O 


Platform 
Displ. 
in 
inches 


-O.20 


-0.40 


Theory —— 
_ Experiment ooo 


Run B1lO Unbraced 


CD 
coos = 1,246 


<< 


9O 180 


ano 


Wave Phase Angle at Front Legs in Degrees 


62 


Run B10Z Braced 


Ud) 
ee ee =a 3 O 
Wn srr? 


VCH 4/19/62 


460 








Figure XIIi - Platform Displacements for Wave C 
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Figure XIIj - Platform Displacements for Wave C 
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Figure XIIk - Platform Displacements for Wave C 
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Figure XII1 - Platform Displacements for Wave C 
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Figure XIIm - Platform Displacements for Wave D 
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Figure XIIn - Platform Displacements for Wave D 
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Figure XIIo - Platform Displacements for Wave D 
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Figure XIIp - Platform Displacements for Wave D 


Theory —— 
Experiment ooo 


0,20 Run D1iO Unbraced 
Platform O wD = 
Displ. | 0 0 Dp, 1.599 





0 90 180 2/0 360 
Wave Phase Angle at Front Legs in Degrees 
0.20 o . Run D10Z Braced 
O (Do 
: Wh, = 5/9 
— 
o 
0 
O 
=-0.,20 


VCH 4/19/62 
70 





IV. DISCUSSION OF RESULTS 


The plots of Figure XII indicate that the results 
of this thesis are good, in the opinion of the authors. 
The procedure set forth in the thesis is judged to be 
sound. 

The reader will notice that in all cases the 
theoretical curve of displacement has the same shape 
as the experimental curve (which is not actually faired 
in). Of particular note are the curves of runs AO and 
AOZ (Figure XIIa), for which both the theoretical and 
experimental curves have an extra dip in the positive 
portion of the curve. 

An attempt was made during the experiment to re- 
late some reference time in the movie and on the San- 
born recorder paper. The procedure involved pushing 
@ manual push-button as the second hand of the clock 
(one rotation = 10 seconds) passed some known time. 

Any human error would appear as an error in the experi- 
mentally determined phase angles 9 at the front legs. 
The resulting data showed that the human error was so 


large that the data had to be discarded. Consequently 
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the phase relationships of the experimental points in 
Figure XII were adjusted to give the best fit to the 
theoretical curve. 

if one follows a particular structure con- 
figuration (for example, 6Z) through its four runs 
he will see that the greatest displacements occur for 
tnose waves which give frequency ratios ¢; of about 1.1 
or smaller. (Note that two different displacement scales 
are used for the plots. ) 

Figure XIII contains one point for each run. Each 
point represents the ratio of the maximum positive 
theoretical displacement, plotted as a function of fre- 
quency ratio # for the particular run. The figure is 
intended to indicate the magnitude of the difference be- 
tween the theoretical and experimental thesis results. 
It is seen that the displacement ratio is generally quite 
different from its ideal value of 1.0. Values greater 
than 1.0 indicate a conservative design, or overdesign. 
Values less than 1.0 indicate that the theory is in- 
adequate and that the structure designed according to 
the theoretical procedure would be too weak. It is in- 
ferred in the previous paragraph that the designer is 
interested in waves for which oa, 18 about 1.1 or less. 


Consequently the right-hand portion of Figure XIII is 


of no further concern, aS a practical matter. For runs 


ne 





Figure XIII - Ratio of Maximum Positive 
Platform Displacements, Thgoretical/Experimental, 
vs. the Frequency Ratio “ for Each Run of the 
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in which it, * less than 1.1 multiplication of the 
theoretical displacements by 1.5 would insure that the 
displacement ratios were always 1.0 or greater. 

The theoretical solution for displacement is very 

sensitive to errors in the structural natural fre- 
quency. Computer calculations (not submitted), using 
both the theoretical and experimental values of Ww, 
obtained in this thesis (Appendix C), indicate that 
in some cases the maximum displacements so calculated 
Giffer by a factor of 2. The computer results show 
that, given two frequency ratios with all other variables 
the same, one may not expect that the frequency ratio 
nearer 1.0 will generally give the larger theoretical 
displacement. If a designer is in doubt about the ex- 
act natural frequency of his structure, as will be the 
general case in design work, he should apply a factor 
of ignorance of 3 to his calculated displacements. 
Under certain circumstances this may be extremely con- 
servative (by a factor of about 5), but at other times 
it may be barely adequate. (The designer may wish to 
apply an additional conventional factor of safety.) 

The authors desired that the theoretical calculations 
Should test the method of determining wave forces on the 


structure and the method for predicting the consequent 
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vibrations. The calculation of theoretical oo was 
felt to be of secondary importance in testing this thesis. 
In order to reduce the sources of errors in the Results, 
tnerefore, the theoretical calculations were made using 
the experimental values of WwW, (as determined in air). 
Theoretical values were used for wave frequency W, 
nowever. 

Additional discussion of results is given in 


Appendices G and I. 
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V. CONCLUSIONS 


ale The procedure set forth in this thesis provides 
a sound method of calculating the dynamic displace- 
ment of a fixed offshore structure acted upon by a 


passing train of regular waves. 


2. The dynamic displacement of the platform may be 
of interest per se (p. 2). More generally, however, 
the dynamic displacement is required for tne design 
Stage determination of the vibration stresses through- 
out the structure (pp. 41-He2). 


5. Field test data on Texas Towers No. 2 (pp. 5-6) 
and No. 4 (p. 1) shows that the highest wave encountered 
does not necessarily produce the greatest platform dis- 
placement. Consequently, for design purposes one might 
not design an adequate structure simply by considering 
only the largest wave to be encountered, as has been the 
general past practice. One must consider several waves 
of various lengths (or periods) and steepnesses over the 
range of possible water depths. The critical displace- 
ment-wave combination must be arrived at by an iterative 


process. 
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4, For a particular structure with its particular 
natural frequency the greatest displacements occur 

for those waves which give frequency ratios a, of 
about 1.1 or smaller. For these waves the theoretical 
solution for displacement is generally conservative, 


with some exceptions. 


ay A designer should apply a factor of ignorance 

of 3 to his calculated displacement. Under certain 

circumstances this may be extremely conservative 

(by a factor of about 5), but at other times it may 

be barely adequate. (The designer may wish to apply 


an additional conventional factor of safety. ) 
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VI. RECOMMENDATIONS 


1. This thesis advances a procedure for obtaining a 
"safe" design of a fixed offshore structure. At 
times the procedure may lead to a considerably over- 
designed (excess strength) structure, which is gen- 
erally poor from a cost viewpoint. It is recommended 
that future workers attempt to refine the procedure 
in order to reduce the amount of overdesign that may 


be involved. 


2. it is recommended that some calculations be carried 
out for some real structures in the ocean and tnat the 
results be compared with observed data. From this work 
one might confirm the validity of this procedure for 
design purposes, and a revised factor of ignorance 


mignt be obtained. For example: 


a; Were the results of the design calculations 
for Texas Tower No. 4 adequate? Could one predict tower 
collapse for the wave and wind conditions that existed 
wnen the tower did collapse? How do calculated displace- 


ments compare for some of the observed data available 


man 2 | ? 
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b. Perform a design calculation for the 
Argus Island structure and compare the results with 


those of [19] : 


c. Perform a design calculation for the 
Buzzards Bay Light Station structure and make a 
conclusion regarding the adequacy of structural 


strength. 


3. Broaden the scope of the design data by extending 
the limits of Figures III, IV, and V to include soli- 


tary and breaking wave theories. 


4, Develop an analogous procedure for calculating 
platform translation and rotation for an arbitrary 


angle of wave impingement. 


5. Continue the investigation of vortex shedding on 

a vertical cylinder in a wave train (Appendix I). In 
a@manner similar to that of this thesis calculate a 
Fourier series representation of the vortex shedding 
force and establish a procedure for calculating trans- 


verse platform displacement. 
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VII. APPENDIX 
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1a 


A. Details of Theoretical Procedure 


Spring Constant k 


Apply simple beam theory [29] to one leg 
(Figure XIV) to determine the static deflection Curve 
for a force Py acting at the level of the platform. 
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Figure XIV - Simple 


Beam Representation 
Gor Sone we SE 


For N legs: k = ie 


Rayleigh's Ener Method A lied to the Unbraced Structure 





Apply Rayleigh's Energy Method [ 22,30 | to the unbraced 


fixed offshore structure to determine the undamped natural 
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frequency Wo: Assume simple harmonic motion for the 
N-legged structure. Use the static deflection curve 


for the elastic curve. 


Woe. 2 «1 GI 
Maximum kinetic energy = 2¢ On Bice +N 2e n max 
0 
2 
2 2 
33° 9g? Wn max 4 12 Nw) 
s - 


Tee Be ie 
Maximum potential energy = N DEL = —SlET 
Oo 


Equate the maximum kinetic and potential energies to 


obtain: 





leg NEL 
a er 3eNwh) (47) 


Use (36) and (46) to find: 


m= 2(W + Z2Nw) (41) 
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If the configuration had been such that the legs 
were considered to be pinned to the platform the fraction 
to be used in the second term of (41) would be different 
from that calculated above. A calculation for this new 
value would proceed in a manner exactly analogous to 
that above. One can deduce that in calculating the 


spring constant k he would start with: 


2 P 
Se grt er f - 2) 


The calculation would give a fraction of 33/140 


to be used in the second term of (41). 


Still Water Added Mass for the Unbraced Structure 


Assume that the added mass for any elemental 
length of leg is equal to the mass of the water dis- 
placed by that leg | 33], measured in terms of weight 
Wy per unit length. Calculate the maximum kinetic energy 
K.E. of this water for a condition of free vibration. 
Then paichiate how much of this added mass unit must be 
placed at the level of the platform (Ww. ) to give the same 
maximum K.E. 
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Use this expression to calculate Wy for d = 27 in. 


oo 


and WV A sleriv 
W = 0.03 N wX 


By comparing this quantity with those of (41) one 
sees that the added mass effect is small. Added mass 


will therefore be neglected. 


Still Water Viscous Dampl for tne Unbraced Structure 





A numerical exanple must be used for this cal- 
culation. Use run O-W, which is the free vibration of 
tne unbraced structure in still water with zero pounds 
added to tne platform in the experiment. Let| y2 denote 
the root mean square velocity of the leg maximum velocity 


profile. Assume simple harmonic motion. 


D = 0.505 in. d = 27 in. R= 42 Ain. 
a 
2 31° 2 
[oe - 9(3Z° — & dZ = 0.1292 y_.¥ 
fo. Ymax/Q2 
O 
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This occurs at Z = 17.0 in., which can be found 
when one realizes that maximum velocity at any point Z 
on the leg is related to maximum displacement at that 
point by the constant wW,° Enter Figure VI with 
a= 17.0 in., a/f = 0.406, ana find Fo/P = 0.360. 

For run O-W the average of the observed amplitudes 
for the decaying exponential displacements was 0.66 in. 

The experimental natural aaeerener was W, = 6.13 rad/sec. 
Assume q and Ww, to be numerically equal, which is very 
nearly true at the small values of damping which are of 


concern. With these numbers one can compute an average 


F iz”) = 0.337 ft./sec. 





Imax = n° Vmax = (6-13)( 


Oe eee 


= 
i 8 = (0.359) (0.337) = 0.121 ft./sec. 


(O. 121) (0.505 505 ) 


which gives a Reynold's number jR = = 4ho, 
(ies 10" } 


Enter Figure II to find steady state Cy - 1.2. 


The maximum damping force acting on the platform is 


p& As Vo where Ag = D.d and N = 4. 


This maximum damping force is calculated to be 
0.00232 lb. At the low Reynold's number involved one can 


assume that the product of drag coefficient C,. and velocity 


D 


84 









is constant, which means that the damping force varies 
Girectly as the first power of velocity rather than the 
second power. This assumption means that the maximum 


damping force equalsC 2 so that 


Le = ; 
Ce = N(C,-¥,) £ D-qe = 0.0192. 


From (41) find m= 0.272. From (37) find 


C, = 2mw, = 3.34. One can calculate viscous a = 0.0057; 
c 


Which happens to be exactly the value that was obtained 
experimentally. It is 11% of the observed structural damp- 


ing (Table X). 


Witn expressions for C and oe one can write: 


ia 
=o N(const. )pD.d. (W,.¥,) 
aa, 





a ha 
U 


For the experimental runs all variables in tne right 
hand side of this expression are fixed except Vo and m. 
(wW cancels out.) If a situation is considered in which 
a7. is held constant, and m is varied, it is seen that 


viscous C should decrease as m is increased. This trend 


C 
Cc 
was observed experimentally. 


85 


Influence Fraction > 


Suppose that the force P acting as shown 
(Figure XV) on a two-legged structure causes the plat- 
form deflection Yinae: The problem is to find what Po 


acting at the platform level will cause the same Ye owe 


> 





Figure XV - Portal Representation of Unbraced Structure 





From (46) it can be shown that: 


PX? 
ymax ™ DAEL and Po = 
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From Figure XV the reader can deduce that: 


P 
PLP - = (sum of forces = 0) 
P 2 | 
and Pa = M, yl (sum of moments = 0) 
For0O#A Za M(Z) = M, - P,Z 
tt Pay = 1 
y" = py = p(y - PZ) 
-) 
yt w AZ - PYZ"y 
EL —— 
2 
ae 5 
ya = eee) 
El’ 2 Ss 
For a = Z =f M(Z) = M, - P,Z + P(Z-a) 


" Ss M, - PZ + P(Z-a) 
Y "=f i us 
. 2 
x vA 
¥' = £7 M,Z - —3— + P(5- - a2) ae 


2 
! a _ Pa 
Equate expressions fory' at Z a to find Cz OET 


= By, 
=k. 1. SF gee ee) ee 
ae - B Ge D 2 








ee 
Equate expressions for y at Z = a to find Cy = — 


Let y' = O for Z-= and use P. =P - 


~o 
at 2 
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Pi 


a 
to find M, Pa(l - By) - = — 


substitute these values for Py and M, into the ex 


pression for y and evaluate Vier aie at Z = ff co find: 


Equate this to the previous value for y from 
Ms ae 


(46) and solve for pe: 
Fo. Bat _ 2a” 
rm “hie p> 


This expression is plotted in Figure VI. 


Fourier Series Representation of Force Acting 
on Platform 


The following two expressions have been developed 


theoretically to represent the force acting on the plat- 


form: 
F(t) = F,sin°wt + Frcoswt (0° wt 4 180°) (55) 
F(t) = -F,sin-wt + F,coswt (180° wt & 360°) (56) 


Many textbooks, such as [22,30], discuss Fourter 


series. 
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let F(t) = a, sinwt + a,sin2wt +... + asin nwt 





2 
+b + b,coswt + b cos2Wwt +... +b cos nwt 
Oo at 2 n 
(O° 2 wit 2 360g) 
a6 20 /w 
where: a= — F(t) | sin not | at 
nS or 
O 
20 /u) 
_ WY 
b= = F(t) | cos nwt | at 
O 
Peary 
58) 
> F(t Jat 
O 
One can deduce that O = a, = ay, a ee 
O = = 5 = oleae 
and by 28 Lr 
ays T/u | a 2 /w 
ere. y “2 g 25 
ay = Sin“wWt dt - — Seoul t edi aq (Fyt i 
Wy 
ow 
sf Joy 
wk WF 
a, = — | sin?w toinseoe] dt - — jsin?wt}/sinswt) 
O w/w 
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al 
25 = tose (Fy + Fo) 


-4 
7 Iso, On (Fy + F,) 


ny 
hh 


-4 
9 6035 (Fy + Fy ) 


Why 2m /uy 
UF, 2 WF, 2 
b, “wiry | sin we || coszwt |at - az | sin wt |]cosaut |at 
uv Oo w/a 
= wa o F,,) 
Dy = O = be = De = 
Thy Om /y) 
OF WF FF, - Ff 
bo = sin*wt dt - 8 sin*wt at = +, 8 
O VU Thy 
nee 
", F(t) = a Sq (Fy + & sin wt 
= tp + F,)sin3Wt + F, coswt 
jor & 2 B, 
F. - F) 
= ee 2"cos2Wt + .... (57) 
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B. SAMPLE THEORETICAL CALCULATIONS 


The calculations performed to determine theoretical 


displacements caused by the B wave are given below. 


Given: 
pile diameter D = 0.5 in. (0.0417 ft.) 
water depth d = 27 in. (2.25 ft.) 
wave height H = 0.313 ft. 
wave length L = 10.86 ft. 
therefore a/L = 0.207 and H/L = 0.0289 


Sample Calculation: 


1. Hydrodynamics 


a.) Calculate maximum drag forces and centroids on 


a single leg. 


oo = 90 <= 


180° 


drag force = Fy RET ams 


Fp = 6.58 Cpa 


where / = 62.4 1bs/ft? 
(12) 


from computer program results (Appendix D) or by in- 


terpolation of Figure IV, A = 0.00421 


on 







SS «<a « Geass 
os © oe, 6 i q 











er iy wae 1¢- hat 





IR = / Sing (18) 


for water temperature = 61°F, y=1.21x 10°? ft.“/sec. 


ee se }> + 39 sinh (25) 

d ~ 2a BS I7H (sinn J) | (0) 
2rd _ 

J = NE = 1.30 


Sinn (2J) = 6.69 , #£=sinh J = 1.698 


0. 0.31 2 + 3(3.14) (0.0289) (6.69 
Te = 26.14) (0.0289 16-69) ] = 0.0715 


8 (1.698) 


i + fe also could be read directly from the com- 


puter program results (Appendix D), 


1+.) = 1.0715 


IR = —O20427 22.2 (age) (00a apni 


ier x 107 1 


from Figure II for /R = 1804, Cy = 0.96 
therefore Fp = 6.58(0.96) (0.00421) = 0.0266 lbs. 


drag force centroid = §'! 


ve 





From computer program results (Appendix D) or by 
ue 
interpolation of Figure Il, == 0.685 


S' = 0.685(2.25) = 1.542 ft. 


ieo° = @ = 360° 


a 


drag force = (F,), es — 0 trough = ie ) 
(Fy )t =~0.58 oe trough = 


from Figure V (or Appendix D), B = 0.00276 





IR >| (19) 
trough a Ua 5 
-H 
EO” 7» Ze -4 (31) 





- H 
or from Appendix D , it + ie = 0.9323 


oy 0.0417 32.2(2 25) (0 00276 } 
/R 5 yf 22:2(2,25) (9.00876) : : = 1591 
trough — 1.21 x 107° 0.9325 
from Figure II for Ry pough = 1591, Cy trough * 0.925 


2 





therefore (F = -6.58(0.925) (0.00276) = -0.01702 lbs. 


pt 


drag force centroid = (S'), 
6") 
from Figure III (or Appendix D), tae 105559 
a 


(§'), = 0.559(2.25) = 1.258 ft. 
b.) Calculate maximum inertia force and centroid on a 
Single leg. 


2 . 
F, ie 7 H tanh J (25) 


for J = 1.30 : tanh J = 0.862 


_ 62.4(0.0417)© (3.14) (0.313) (0.862) 
. ee 


F = 0.02295 


al 


= _ oD? [1 +5 sinn J - cosh J | 
i 8 Fr cosh J 


sinh J = 1.698 F cosh J = 1.971 


5 ey 0.0417 ° 10).86))) (One 7, 1 + 1.3(1.698) - 1.972 | 
Je 0.02295) (1.971) 


S_ = 
; = 1.26 ft. 
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therefore (F = -6.58(0.925)(0.00276) = -0.01702 lbs. 


pt 


drag force centroid = ($'), 
(S') 
from Figure III (or Appendix D), t = 0.559 
d 


($'), = 0.559(2.25) = 1.258 ft. 


b.) Calculate maximum inertia force and centroid on a 


single leg. 


2 : 
F, = = w H tanh J (25) 


for J = 1.30 , tanh J = 0.862 


_ 62.4(0.0417)* (3.14) (0.313) (0.862) 
oe) ae 


F = 0.02295 


i 


Y D° LA [ 1 +5 sinh J - cosh J | 


Sr * 3 F, cosh J 


sinh J = 1.698 : cosh J = 1.971 


5 . 62.4(0.0417)° (10.86) (0.32 [3 + 1.302.698) - 3.972 | 
I 0),02295 ) (1.9745) 


S_ ss 
1 7 1-26 ft. 
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ec.) Calculate wave period. 


. d_ tanh J$1 ef 2(cosh(2J))*+2cosh (25) +5 
T = C where C = ga tanh J} om 18 — on a ) h 


(3) 


C _ 222 2.25)(0.862 f2+(3.14(0,0069))" 216.11)" 22(6.77)45) 
8 (1.698) 


Cc = 6.98 ft./sec. 


Therefore T= Las = 1.558 sec. 


Doe = = 4,03 rad./sec. 


C is determined internally in the computer program 
(Appendix D) but is not one of the outputs. A simple 
addition to the "WRITE OUTPUT," “DIMENSION,” and "FORMAT" 


statements would include C as an output. 


Summary of steps a,b, and ec: 


Fy = 0.0266 lbs. S' = 1,542 ft. 
(Fp) = -0.01702 (OM mel,250 ft. 
F_ = 0.02295 lbs. S.-= 1.26 ft. 


I 
WwW = 4,03 rad./sec. 


ig 
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Enter these figures in the upper left-hand corner 
of the tabular ealcullgcaem ameen. Figure XVI. (Actually 
in this sample calculation the numbers do not ieee ex- 
actly with tnose of Figure XVI, which were determined at 
another time, but the slight differences are negligible 


for purposes of illustration. ) 


2. Vibrations 


Consider run B6. This is the unbraced model 
with 6 lb. of added weight on the platform, in wave B. 

The reader should trace the steps on Figure XVI 
as he reads the following explanatory notes. 

Values of WW? —, and k are taken from 
Appendix F. Otherwise Biey would be determined as ex- 
plained in the Procedure, Section II. | 

Values of the influence fraction =o are read 
from Figure VI. 

The values of QO ‘at the front legs are simply 
reference values. Values of platform longitudinal 
displacement are calculated for these twelve reference 
values of 9, and a curve is faired through the points 
(as was done in the preparation of Figure XII). 

Below the values of 9 the corresponding re- 


ference angles for each of the significant terms in 


96 












—_— ——- =-_-— = ——- — ame 


_— | *—— > a cae § 
° = | eee eee 


eo SS _ @& ——_ es =f oo = i) 
= ss ' a § GHP) equ Gage 
=) a \oe 


- = =a at oe & > <-> 


| 


—_— = — ee ee 
aa eet eS aw 


—s 2° — wey Ss —, 


a ee 








eet _[ eset [aco [bo2 [ce [Lvs Tes" [Sve [@00] [sw] Coven) anavanrusa varaeaiono7 ware 























}o0’- j}o0°-| oO | /00' 
Lai'- v0C' |922' 1931! | 001' |9I0'- 
l0g'— loo'-|100'-|100' |100'| Joo" 
9g1'- ‘| 280|-400'-|So/'—|b91'~ 88 |'—|| 


O | (00) (00-| O | /00'| 100'| © | joo'-| 1o0'— O 
$ll-|¥22-| vIZ-| Iel'-| Leo-| 820°) +zf'| pez] vlZ) ont: 
0 | 100°) O | lo0'-| O | 100’! Oo | Io0'-} @ | Joo’ 
IZI—|¢¢o'-| 90°} eel) S8t'| LL] ta) sgo'l po0'-] aer- 


Bee [ 28 oe | 900" se"| a0 Lee [oo |e" | a [0 0 30 
aS SS eee 


ABIL ILI |L |L0E|2¢2|.281|.221|29 |2 |zo¢ _/me) 
92E| 962 at IOCILI1 II oe 295 |.9Z 
of €2) Jol lis | IE) EZ) Illi | Ze |.1€Z 
aid cf 


Xo 


| | 





FOZIILI|I¢! 1.911 |.98 |.9S 9 4. Ibe 
SSI|Sb ISé€ |s¢eeled2|.sl2l|.acl |b |oe 
OE |.01€ | 082)052].022|.0b! | 091] .0€1|.00! |.02 Yen S*228/av4 SO p= 1m 
eh |.€ jecc|.esl|eb |e |ed2| £9 |b |.¢ 8800 OF? asta go'y = 
OC?| O41] 091 |0¢! |00!|.02 | oF |.0!1 |.O+€ 
Oise Pitt ee 99 NOY wos LNAIWZBI071d51¢ 
ee ee ee ee ee TYNIdnLt9NO7T Wear Td 


Se ne nnn en NE vn 


O42} 91d 
OL 1.09 |.0€ |.09¢| OFE|.90£f| .0472|.0b2| .0!1Z).091|.051 
0O9E 1062) .08/|.96 |.09¢|.02Z| 08] |.0b | .09¢|.922|.081 
CIS | OEE | OOF |.0LC | .ovCc| O12] .08f |.051|.92! | .0b | 09 


09F 88 |ooe OL2| 02 212 0a r 
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the Fourier series (57) are listed. It is apparent that 
a cosine varying force term leads the sine varying force 
term of the same frequency by a 90° time phase angle. 

In order to fit the cosine terms to the classical dif- 
ferential equation of forces (32) 90° is added to each 
of the cosine term reference angles, after which the 
term is treated as a sine function. 

In the middie portion of the table the actual 
Va.ues of the angles for each of the terms are listed. 
For the front legs the angle @ is determined from (40) 
Seetable Il fer each velue of RW, where n is the co- 
efficient of wt. The angle @ ie subtracted from the 
reference angle for each case, because displacement 
lags force (33). At the rear legs the force lags the 
force at the front legs by the lag angle nZ,, (58). 
For each specific case the angle at the rear legs is 
found by subtracting nar, from the corresponding angle 
at the front legs. 

In the lower portion of the table the values of 
e (33) are calculated for each term of the Fourier 
series for displacement. The expression for the force 
coefficients F are given in (57). Values of the mag- 
nification factor em are determined from (38) 


ASTATIC 
or Table I. 
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The reader is cautioned that linear interpolation 
of Tables I or II may be poor in the immediate vicinity 
of resonance iF 0,). 

In each of the twelve columns the value of Xo is 
multiplied by the sine of the appropriate angle from 
above, and the result is recorded. These results are 
added in each column to give values of platform longi- 
tudinal displacement x. 

For this sample calculation the ratio of the maxi- 


mum positive theoretical displacement to the maximum 


positive experimental displacement is 1.17. 
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C. Summary of Theoretical Calculations 


Forces and Force Centroids 


Table IV describes the 4 waves used in this thesis. 


The values H and L are measured, while a ito 


is read 


directly from the computer output in Appendix D (or from 


(6)). For all runs d = 2.25 ft. 


Frequency (radians/ 


second) is calculated from (8) by use of the relation 


CO = 29rC/L. 


TABLE IV - Wave Cnaracteristics 


sf 49. -H 
Wave H(ft.) L(ft.) @/L H/L J=2nd/L TS ee 


A 0.240 
B OVS13 
© Creon 
D 0.404 


1Dm2o 
10.86 
See 


(og 


0.1476 0.01575 0.928 
0.207 0.0289 1.300 
0.255 0.0409 1.603 
0.307 0.0552 1.930 


1.0557 
Odes 
Ooi 
1.0908 


0.949 
0.932 
Oyoes 
Q.911 


Table V lists parameters used for computing drag 


forces. The multipliers A and B are read directly from 


the computer output (or from Figures IV and V, re- 


spectively), IR and Rerougn are computed from (18) and 


(19), respectively, while C, and C 


Weeare Il. 


100 


D trough 


are read from 


CU 


Blue) 
4.00 
4.60 


5.13 
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TABLE V 


- Various Coefficients for Drag Force Determination 


Wave A B IR Ch trough Cy trough 
A 0.00272 0.00191 1410 0.846 1188 0.846 
B 0.00421 0.00276 1804 0.925 1591 0.885 
C 0.00538 0,005 2070 0.925 1720 0.925 
D 0.00674 0.00342 2290 0.962 1795 0.925 

Table VI lists the theoretical forces and force 
centroids on a single leg of the structure. Fp and (FL), 


are computed from (12) and (13), respectively, while Fy 


is computed from (25). The centroids S' and’ (S'), are 
read directly from the computer output (or from Figure III), 


and Sr 


measured from the bottom of the tank. 


is computed from (28). Centroidal distances are 


TABLE VI - Theoretical Forces and Force Centroids on a 


Single Leg 
Wave *D(ibs) Si (rt) (Fy), (1bs.) (S'), (ft) F,(1bs) °r (ft) 
A 0.01755 1.375 -0.01l242 1.180 0.01485 1.200 
B 0.0266 1.542 -0.01702 1.258 0.02295 1.260 
c 0.0337 1.680 -0.0203 1,550 0.0283 ibys 
D 0.0425 1.813 -0.0216 0.0329 1.380 


1.408 
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Structure Natural Frequencies 


Theoretical natural frequencies of the structure 
are computed by (47) and (61), using numbers from 
Appendix E. In Table VII run 3, for example, designates 
the unbraced model with 3 pounds of added weight on the, 
platform. The letter Z designates the braced model. 
Experimental values of Ww (Appendix F) are also given 


in Table VII for convenience of comparison. The units 


of all frequencies are radians/second. 


TABLE VII - Structure Natural Frequencies 


Run Wn theoretical “n exp. in air Wn exp. in water 
0 5.25 6.01 ~ (es 
3 453: 4.67 4.79 
6” 4.05 4.00 4.02 
10 3.59 Seek One T 
OZ 5.20 5.97 6.09 
52 4.50 4.95 5.00 
6Z 4.02 ae 41 
10Z pea Sees Bel 
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D. Computer Programs 


Two separate IBM 7090 FORTRAN computer programs were 
written for use with this tnesis. Part 1, "MAX DRAG 
FORCE MULTIPLIER AND DRAG CENTROID,"” was used both to 
develop the curves families of Figures III, IV, and V, 
and to supply exact values for force and centroid de- 
termination for the 4 waves studied. 

Part 2, "TOWER DISPLACEMENTS," was used to supply 
the theoretical platform displacements for the 32 ex- 
perimental runs of this thesis. The output from this 
program is snown plotted as Figures XIiIa-p. The output 
from Part 1, although comprising data for 76 waves of 
differing characteristics, is included in this Appendix 
only for those waves wnich were the object of experi- 
mental study in this thesis. (The remaining output data 
has been submitted to the Thesis Supervisor under 
separate cover.) 

Because the designer may want to make use of the 
programs, a brief description of their characteristics 


is included. 


Part l 


This program was developed to facilitate solutions 


of drag forces and drag force centroids for crest and 
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trougn regions. It is applicable only for those wave 


parameters for which Stokes! 4rd approximation applies. 


Its primary use is for development of curve families 


on (S! )t 


such that A, B, 5 and a: aa in tne formulas listed 


below may be read directly by the designer if no com- 


puter is available. 


» = & da* “Da 
D 2 
_ypdac 
(F ) = D trough B 
DG 2 
1+ o/a 
ci - S' — max a(S'/d) 
cy a ed 
O 
A 
1+ o-H 
. 2 
was S' “trough d(S'/d) 
ok a ged 
O 
B 
1+ o/a 
2 
where A = “max.a(S'/d) 
exe | 
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(12) 


(13) 


(26 ) 


(27 ) 


(14) 






14+Vo-H aa H 


ef Yerous d(S'/a) (15) 


Integration of (26,27,14, and 15) is performed 
by the trapezoidal rule with 21 unequally special 
stations. 

The curve families developed (Figures III, IV, 
V) show lines of constant H/L on graphs whose co- 
ordinates are the desired values vs. d/L. 

A secondary function of the program is to supply 
horizontal-component velocity profiles at wave angles 

= 90° and 270°. Two uses of these profiles are to 

give specific indication of the limit of applicability 
of tne authors! design procedure and to supply required 
information for an analysis of transverse vortex 
shedding forces (Appendix I). 

Supplementary information which may be read directly 


from the output includes: 


a+ No a + o-H U nex UE rough 
3 3 3 
a a Vea Vea 


LOD 


: 
| 








Wave celerity C, although determined internally as 
Fae not included in the output. However, the 
ll inclusion of this parameter (called "T" in the 
program) in the WRITE OUTPUT TAPE, DIMENSION, and out- 
put FORMAT statements would cause it to be written 
output. 


Variables to be read into the program on data 


cards are defined as follows: 
X= di 
Y = L/H 


The above data for each run is put on the same 
card in accordance with FORMAT statement 10. Variables 


written as output are: 


X = G/H 
Y = L/H 


qa 


PC = (S'/d)n n= 0, 20 
u 
uc = (=), n= 0, 20 
ga 
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a 


u 
ucs@ = (3), nia 'O, 20 





S! 
pre (—SECush) n= 0, 20 


u 
Wie ( proueny 
Yeh 
7 2 
UTSQ = (EEE) n=0, 20 


Values of PC, UC, UCSQ, PT, UT, and UTSQ are given 
in incremental steps of S/d (not S'/d) of 0.05 from 
the bottom to the still water level. At n= 20 it is 


seen that the following result: 


a+“7o 
oe co 


U 
Uc = vex 


y gd é 


U 
UCSQ = = 


- H 
aU ae ) 
PT = a 





U 
UT = trough 


Y gd 
2 


Ue rough 


UTSQ = zd 


The output is written in accordance with FORMAT state- 
ments 90 and 100, and is included in this appendix for 


the A, B, C, and D waves only. 
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If other than 36 input data cards are to be used, 
the program must be run separately for each set of 
36, with blank cards being used to fill out any set 
not comprising 36 data cards. A card with ene 
punched in the x columns should be used to separate 
the data from the blank cards. 

Some of the interesting aspects of the output used 
to develop Figures III, IV, and V are discussed in 


Appendix G. 
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* M1933-1598 , FMS , DEBUG ,2,4,1000,0 


J XEQ 
* eST 


W.C. NOLAN. AND V. C. HONSINGER 


WAVE INDUCED VIBRATIONS IN FIXED OFFSHORE STRUCTURES 
PART 1, MAX DRAG FORCE MULTIPLIER AND DRAG CENTROID 


C262 CQ 


DIMENSION X(36),Y¥(36),A(36),B(36), (21,36), 08 er ae »PC(2lIe6)e 
IPT (21 , 36) ,UC (21,3 ) ur (2 ,30),UCSQ (21, 23 ), UTSQ (21 
READ INPUT TAPE 4,10, (X(I), Y(I), Tl, 
10 FORMAT (2F10.4) 
DO 60 I=l1 ,36 
TF (X(Z)) 
IF (X(I)) 70,70,20 
20 Z=6, 2832*x (L)/¥(I 

D=EXPF (Z )-EXPF (-Z 
E=TANHF (Z ) 
FeEXPI' (2. *Z )+EXPF (-2.4*Z) 
T=SQRTF ((E/Z )*(1.4+(9.8697/¥ (L)**2. )*(F**2.4+2.*F+10. )/D**4, )) 
R=0 
DO 50 J=#1,21 
GP=EXPF Ps ae 
GM=EXPF (Z*R )-EXPF (-Z*R 
HP=EXPF (2, *Z*R)+EXPF (-2. *Z*R 
HM=EXPF (2. *Z*R )-EXPF (-2.*Z*R 
PP=R+9. 4248 *HM/(2.*X(I)*Y (I) *D*¥*4, ) 
UU= (9, 8697*T/ ( Y (L)#**2, )*D**2,))*(-2.46. *HP/D**2. +HP ) 
PC (J, I )=PP4+GM/(2.*X(I)*D) 
PT (J,1I )=PP-GM/(2.X(1I)*D) 
UC (J, 1 )=UU+3. a 

| UCSQ (J I )=UC (J,1)**2. 

ire (sy, TauveS. viaoensée/(x(t)*D) 
UTSQ (J, I )=UT(J,1)** 
IF (J- 7) 30,30, is 








30 A(I)=0 
B(I )x0 
WC=0 
WT Q 
GOTO 50 
LO - en au bee ee EA ee (2. 
B(I)=B(I )+(UTSQ(J,I )+UTSQ(J-1,1))*(PT(J,1)-PT(J-1,1))/2. 
eat A RR Meese ta le ae 
Wl=WT +(UTSQ (J, I )+UTSQ(J-1,1))*(PT(J,1I)**2.-PT(J-1,1)**2. )/4. 
50 R=R+.05 
Peas (view /d (x) 


60 CENTT (I )#WT/B(I) 

70 DO 80 Mel ,N 
WRITE OUTPUT TAPE 2,90,X(M), 
ee M),UC(1,M),UCSQ(1,M),PT 
sg E gurpyt TAPE 2, 109, (PC ( 
1) UEee (J 


Jee ol 
80 CONTINUE 
90 FORMAT Sn cee 
100 FORMAT (56X ,6F9.6) 
CALL EXIT 
Boden ,0,0,0,0,0,0,0,0,0,0,0,0,0) 109 


(M),A(M),CENTC (M), BM ,CENTT (M), 
(1,M),UT(1,M),UTSQ(1,M 
js UC (J,M), hit bS M), PT(J,M),UT(J,M 











— Gee 6 
hw 
=~ - 
; «= 
=> -_ => _- =& Le 
~ 
‘5 
_ —=_= oe 
= ———— 


ae] 
La thli aM ‘/ y ® } 








, anes F 





iis | 
| + | B St 

| Pees a eee 
ime 00000 0.002722 gre icelazonsz 0.524623 


t 








0. 
0.052390 








gd 








cae seeee 
trough | 
e eo ea 6 
0-042524/0.001808!0. 
0.042578! 0.001813'0.04775 













“+ rough , u 


cea | aa 


0.040005 ‘0.001600 
0.040039 0.001603 














0.104785 |0-042740 0.001827 | 0.095500-0-040144 0.001612 

0.157191 |0.043011)0.001850 0.143234%-0.040319 0.001626 

0.209614 | 0.043391 /0.0018483 0.190965-0.040563 0.001645 

0.262060 |0.043882 0.001926 0.23867 3-0-040877 0.001671 

0.314534 |0.044486 0.001979 0.286 359-0.041262 0.001703 

0.367043 |0.045204 0.002043 0.334019-0.041717 0.001740 

0.419591 |0.046040 0.002120 ' 0.381648-0.042242 0.001784 

ee 0.472186 |0.046994 0.002208 0.429241-0.042837 0.001835 
0.524833 |0.048072 0.002311 0.476794-0.043503 0.001893 

0.577539 |0.049277 0.002428 0.524303-0.044240 0.001957 

0.630310 !0.050612 | 0.002562 0.571762-0.045047 |0.002029 

0.683153 /0.052082/0.002713 0-619166-0.045926 0.002109 

0.736075 |0-053692 0.002883 0.666512-0.046875 0.002197 

0.789083 10.055449 0.003075 0.713795-0.047895 0.002294 

0.842185 /0.057358 0.003290 0. 761004-0.048985 0.002400 

0.895387 |0.059425 | 0.003531 0-808148-0.050145 0.002515 

0.948699 [0.061659 | 0.003802! 0.855209-0.051375 0.002639 

1.002129 0.064067 0.004105 | 0.90218 -0.052675 0.002775 

114055685 [0.066658 | 0.004443! 0.949075-0.054043 ,0.002921 

Be eo c00 [000828 lovéesean[o.cozrse 0.559225!0. 0-044294/0.001962|0. 0.043083 |0.001856 
0.052704 10.044401 |] 0.001971] 0.047370-0.043160 |0.001 863 

0.105426 10.044724 | 0.002000} 0.09473 0.043393 |0.001 883 

| 0+ 198160 | 0.045263 | 0.002049 0.142064-0.043780 |0.001917 

10.210935/0.046024 0.002118) 0. 18937 -0.044323!0.001965 

0.26375910.047009 | 0.002210! 0.23663 7-0.045024 }0.002027 

0.316644 | 0.048225 |0.002326| 0.283 844-0.045882 |0.002105 

0. 369602 |0.049680| 0.002468! 0.3 30990-0.046900 !0.002200 

a |0-422646 0.051382 |0.002640! 0.37805 -0.048078 |0.002312 
|0-475791 0.053342 |0.002845 | 0.425034-0.049419 0.002442 

, 0+ 529050 | 0.055572 | 0.003088 | 0.471915-0.050925 |0.002593 

' 0.582438 |0.058087 | 0.003374] 0.51868 I-0.052596 | 0.002766 

(0.635970 10.060901 0.003 709 | 0.565 322-0.054435 | 0.002963 

0-68 9663 |0.064033] 0.004100! 0.611824-0.056444 | 0.003186 

0. 743532 |0.067505 | 0.004557 | 0.6581759-0.058623 | 0.003437 

0.797597 (0.071337 | 0.005089: 0. 704 360-0.060974 '0.003718 

0.851875 10.075557} 0.005709: 0.750 364-0.063497 | 0.004032 

0.906386 {0.080192 | 0.006431 | 0.796174-0.066192 | 0.004381 

0.961152 }0.085275 | 0.007272 0.441 769-0.069059 | 0.004 769 

1.016196 |0.090841 | 0.008252 | 0. 887139-0.072095 | 0.005198 

1.071540 |0-096930] 0.009395! 0.932264-0.075296| 0.005670 

62230000 124.450000 10.005 381 |0.746012 {0.003177 [0.5907421 0. 0.041733!0.001742/0. 0.041139!0.001692 
| | | | | | 0.052724 10.041884| 0.001754 | 0.04731 9-0.041254/ 0.001 702 
0.105465 |0.042338| 0.001793 / 0.09461 ¥-0.041598! 0.001730 

0.158243 |0.043100| 0.001858) 0. 141878-0.042173{ 0.001779 

0.211074 |0.044176, 0.001951! 0.18909l-0.042981 | 0.001847 

0.263979 0.045575 1 0.002077 | 0.2 36236-0.044024 0.001938 

0.316974 |0.047310' 0.002238 | 0.283295-0.045306 | 0.002053 


0.370081 0.049396 0.002440! 0.330250-0.046830!0.002193 

















mee Beau 0.002689 0-377081-0.048601 | 0.002362 
~ oaee 0.476713 102054702] 0.002992] 0.423 76 4-0.050623! 0.002562 
0.530282 |0.057970] 0.003361 | 0. 4 7029 4-0.052902 | 0.002799 

0.584050 | 0.061688] 0.003805 | 0-516633-0-055442 | 0.00307% 

0.638043 | 0.065890] 0.004341 | 0.56276 I-0.058249/| 0.003393 

0.692287 |0-.070616]| 0.004987 | 0.608664-0.061327, 0.003761 

0.746812 102075914 | 0.005763 | 0.654314-0.064680! 0.004183 

0.801648 | 0.081835 |0.006697|0.699675-0.068310, 0.004666 

0.856828|0.088440| 0.007822 | 0. 744724-0.072220 0.005216 

0.912388 |0.095798 | 0.009177 | 0.78944 3-0.076407 0.005838 

0.968365]0-103985 | 0.010813] 0.833784-0.080869, 0.006540 

1.024800 | 0.113090 |0-.01278910.877724-0.085597' 02007327 

1-081740 {0.123214 0.015182] 0.92122 7-0.090581 | 0.008205 

nei ge |0-006742 (0-806107 [0-003618 lo. 626576 o 0.036887 |0.001361/0. 0.036639 0.001342 
0.05257810-.037079 | 0.001375 | 0-047440-0.036791/ 0.001354 

0.105180! 0.037655 0.001418 | 0.094854-0.037246| 0.001387 

0.157831 10.03862410.001492 | 0.142224-0.038008| 0.001445 

0.210556 | 0.039998 [0.001600 | 0. 18952 L+0.039081| 0.001527 

0. 263380!0-041792 10.001 747 | 0.23672 1+-0.040470| 0.001638 

0. 316330} 0.044031 |0.001939|0.283798-0.042182)| 0.001779 

0. 369435 10.046741 {0.002185 | 0.330726 0.044227/ 0.001956 

57 024227241 0.049958 0.002496 10.377477-0.046614| 0.002173 
= 0.476227 10053723 [0.002886 |0.424020-0-049354| 0.002436 
0.529978] 0.058086 |0.003374| 0-4 703250052458! 0.002752 

0.584013] 0.063106 [0.003982 10.516 357+ 0.055937| 0.003129 

0.638371] 0.068849 {0.004740 | 0. 562080+0. 059801 | 0.003576 

0.69309310.075398 [0.005685 | 0.60745 6-0. 064062! 0.004104 

0.748225! 00082844 10.006863 10.652442+-0.068725 | 0.004723 

0-80381610.091297 [0.008335 [0.696994+0.07 3796 | 0.005446 

0.859920| 0. 100882 ]0.010177) 0+ 74106 3-0.079275; 0.006285 

10.916597|0.111747 [0.012487 |0.784596-0.085155! 0.007251 

0.973913]0.124064 0.015392 | 0.827536--0.091421; 0.008358 

1.03193810.138031 [0.019053 |0.869820 0.098046 | 0.009613 

1.090752| 0.153882 }0.023680 0.911 380F0. 104989, 0.011023 


110 








PART 2 


This program performs the tabular calculation of 
Figure XVI. It is written for the specific case of 
N = 4 legs and k = 0.625 1b ioe In the expressions 
for 5, T, U, V, and W the constants are recognized 
to be x or x For different values of N or k the 
constants may be changed by substituting new values 
in the program. 


Variables to be read into the program on data 


cards are defined as follows: 


y = run number 
eC) 

A Lor 

B= c/C 

C = Dror (radians ) 


- F 
pp = —-,—*, (1bs.) 


_4 (F, + F,) (ibs.) 
P = oe ( 1 2 
4 (FF, +F,) (1bs.) 
a ise Le 
R = F, (lbs. ) 


The above data for each run is put all on the 


same card in accordance with FORMAT statement 10. 


14D 





Variables written as output are: 


y = run number 


x = platform longitudinal displacement (in. ) 


This output is written on a single line for each run, 
in accordance with FORMAT statement 100. For each in- 
put data card 12 values of platform displacement are 
computed, these values corresponding to the le angular 
wave positions 9 indicated at the top of FPigure XVI. 
If more than 32 runs are to be used, the program 
must be run separately for each set of 32, with blank 
cards being used to fill out any set not comprising 
32 data cards. 
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*M1933-1598 , FMS , DEBUG ,2,4,500,0 


20 
30 


HO 
50 


60 
(0 


1ETA(3,J)-F) + W*SINF (THETA (4,J 


XEQ 
LIST 


W.C.NOLAN AND V.C.HONSINGER 


WAVE INDUCED VIBRATIONS IN FIXED OFFSHORE STRUCTURES 
PART 2, TOWER DISPLACEMENTS 


DIMENS LON ¥(32),A(32) ,B(32),C (32) ,DD (32) ,P(32) (32) R(32) THETA (4 


1,12) ,X (32,12) 


READ INPUT TAPE 4,10, (Y(I),A(I),B(I),C(I),DD(I),P(I),Q(I),R(I),I=1 


se): 
10 FORMAT (IF 3.0,7F8.5) 


DO 90 I=] ,32 

D=2.*C (I 

E=3.*C (I 

FFe2.*A(I)*B(I)/(1.-(A(1))**2. ) 

IF (FF )_20,30,30 . 

PR=PF+3.14159 

CONTINUE 

F=ATANF (FF) 

GGw ,.*A(I)*B(I)/(1.-4. *(A(L) )**2. ) 

IF (GG) 40,50,50 

GG=eGG + 3.14159 

CONTINUE 

GeATANF (GG) 

Bee st on T)/(1.-9. (R(T) eee) 

IF (HH) 60,70,70 A 

HHsHH + 3.14159 

CONTINUE 

H=ATANF (HH) 

AAzw1 ./SQRTF 

BB=1 . /SQRTF 

CC=1./SQRTF 

S=6 .4*DD (I ) 

te 2 

Um . 2*Q aot 

Ves. 2*R(I)*AA 

Wa -%,2*DD (I )*BB 

DO 80 J#1,12 

Ze J 

THETA (1, J )#0.5230*Z 

THETA (2, J )#1.5708*Z 

THETA (3,J)=1.5708 + 0.5236*Z 

THETA (4,J7)#1.5708 + 1,0472*Z 

X(I,J)aS ai a a iH + U*SINF (THETA (2,J)-H) + V*SINF (TH 
-G) + T*SINF (THETA (1, J sie t + U 


eereartea tec tec sie 


T)*B(L))**2, 
A ¥HD, 


*A(I me 
6.*A (I )*B(I 











<a 


HED, 


1.-9.* L ERO, Oe 


O*SINF (THETA (2, J)-H-E) + V*SINF (THETA (3,J3)-F-C(IL)} + W*SINF (THETA (4 
3,J)-G-D) | a 
80 CONTINUE 
90 CONTINUE 
100 FORMAT (1F6.0,12F7.3) 


110 


DO 110 I=#1,32 

WRITE OUTPUT TAPE 2,100,Y(L), (X(I,J),J=1,12) 

CONTINUE 

CALL EXIT i 
END (1,1,0,0,0,0,0,0,0,0,0,0,0,0,0) L135 





E. Sample Calculations for Experimental Data 


B wave 


Wave period T: 
Sanborn recorder speed = 25 mm/sec. 
Wave length as measured on recorder = 40 mm. 


Therefore T = 40/25 = 1.60 sec. 


Wave circular frequency w: 
ow or 
oo: =“ 740" 3.93 radians/sec. 


Wave height H is read from the experimental wave profile 
(Figure XVIIb). 


Wave forces: 
Wave drag and inertia forces are read directly from 


the experimental force record (Figure XVIIb). 


Fy (9 = 90°) = O.O44 lbs. 
(Fp )+ (0 = 270°) = -0.039 lbs. 


F (9 = 180°) = 0.017 lbs. 


ii 


Properties of the Structure Model 
Consider run 6 and calculate the structural 


damping by (42) and (44). 
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everege Si = G2 + 2 + 6b 4 2 AB) wai 


d= log, (1.47) = 0.385 (42) 
C _ 0.38 
= es 0.0613 (44 ) 


Cc 


Damped natural frequency: 


q = Stes (an) = 3.99 radians/second 


From (45): 


CU eee 4.00 radians/second 


n 
fi - (0.0613)* 


The volume of plastic in the unbraced model (in- 


cluding base) is: 


. 21 a sal 1 
asd)" ; £(1.)°|| 0-619] + ye (0.505)" (42) + (30g - q) (22) (5) 


#199 + 33.6 + 329 = 561.6 tm 


For one leg: 


wR = (4) (22:9e) (20h) 10.3H0 ad... 


LL 





For the platform: 


The bolt weighs 0.200 lb. 


W = (529-5) (228) ‘+0 -2007008.225. 2b. 0 


Calculate the following non-dimensional ratios. 


From (59): 


_ €4)(0.340) + 0.4 
ty 7 ay tossHoy 7 2-29 


From (63): 
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F. summary of Calculations for Experimental Data 


Wave Characteristics 


Wave profiles for waves A-D are:shown in Figure XVII 
a-d, respectively. The crest is in all cases the 8th 
wave crest to pass the cylinder. Wave phase angle 0 


is set equal to 90° at the crest. 


TABLE VIII - Experimental Wave Characteristics 
st (oe) 


Wave H(ft.) L(ft.) d(ft.) (séc. ) (radians/sec ) 


A 0.240 15.23 2.25 2.00 3.14 
B 0.313 10.86 2.25 1.60 5598 
¢ 0.361 8.82 2.25 1.40 44g 
D 0.404 135 2.25 1.26 4.99 


Longitudinal Wave Forces 


Experimental and theoretical longitudinal forces 


associated with waves A-D are shown in Figure XVIIa-d, 


respectively. 


F. at 0 = 90° 


D 
(Fp), at 9 = 270° 


= ° 
Fy at 9 0 
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The forces are read as follows: 
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Figure XVIIla also shows experimental transverse 
forces caused by vortex shedding in the A wave. 
The experimental longitudinal forces read from 


Figure XVII are given in Table IX. 


tt 
TABLE IX - Experimental Wave Forces, 5 Cylinder 


wave "p (1ps.) ("o)t (10s. ) Fl (ibs. 
A 0.014 -0.019 0.017 
B 0.044 -0.039 0.023 
C 0.060 -0.055 0.022 
D 0.085 -~0.039 0.032 
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Figure AVIIb 
08 Longitudinal Forces 
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Figure XVIIc 
Longitudinal Forces 


1/2" Cylinder - C Wave 


Z Theoretical 


Oo oO O Experimental 


Wave Profile - C Wave 





0-0-0 Experimental 
WCN 4/13/62 


50 00 10. 2 2298 20 60 0 
421 


i 


Phase Angle in Degrees 





“10 Figure XVIId 
.08 oa ° LONGITUDINAL FORCES — : 
1/2" Cylinder - D Wave ° © 


Theoretical 





06 0 0 O Experimental 


Wave Profile - D Wave 


Wave 
Profile 
in Feet / 





0-0-0 Experimental 


Phase Angle in Degrees WON 4/13/62 


— 


‘) ho 3: 390 160 200 240 280 D0 360 5 S¢ 0 


122 





Properties of the Structure Model 


TABLE X - Experimental Damping and Structure Natural Frequencies 


Run Total Structural Viscous Damped Undamped 
Damping damping damping Natural Natural 
(as fraction C (as fraction Frequency Frequency 
of critical) e of critical) q(rad/sec) 4 (rad/sec ) 
O 0.0501 0.0501 0 6.00 6.01 
5 0.0548 0.0548 O 4.66 4.67 
6 0.0613 0.0613 O 3.99 4,00 
Oo 0.0748 0.0748 O Feo Daek 
O-W 0.0558 CTOpen OulOaws 6.12 O15 
3-W 0.0580 0.0548 0.0032 4.78 4.79 
6-W 0.0655 0.0613 0.0042 4.01 4.02 
10-W 0.07438 opney gts O 4.26 Soar 
OZ 0.0501 0.0501 0 5.96 Son 
py 0.0559 0.0559 O 4, 94 4.95 
6Z erroleney) 0.0613 O 4.09 4.10 
10Z 0.0738 Cm OW ears: O Beet 3-25 
OZ -W 0.0547 O75 C1 0.0046 6.08 6.09 
3Z -W 0.0592 0.0559 OnOO55 4.99 5.00 
6Z -W 0.0645 Opnelon O20052 4.10 ‘Teen 
10Z -W ORO50O 0.0738 ~- Beer 3.44 





Spring constant k (46) was determined experimentally: 


k = 0.625 pounds/inch 
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G. Supplementary Discussion of Results 
1. Theory and Experiment 


Longitudinal Forces 


A glance at the longitudinal force plots of 
Figures XVII a-d shows considerable disagreement be- 
tween theory and experiment. Theoretical forces (drag 
and inertia) are tabulated in Table VI, and experimental 
forces are in Table IX. These tables show close agree- 
ment of inertia forces, with the drag forces differing 
by a factor of 2, in cases. 

Prior users of the force gage have indicated that 
the dynamic response of the gage to exciting forces 
may indicate forces greater than the actual [es]. 
However, the authors conclude that for the frequencies 
associated with the waves used in this thesis any magni- 
fication factor applied to force gage displacements | 
would in all cases be less than 1.1. 

The proposal to use tne rms: velocity in the Reynolds 
numbers for Grag coefficient determination is based upon 
the assumption that Ch and Cy trough are independent of 
depth [27]. For the waves studied the lowest velocity 


encountered in the crest region is (theoretically) at the 
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bottom of the tank for the D wave. From Appendix D 

it is seen that “max = 0.0369 for this situation. 
Therefore, Vea uo 0.315 ft./sec. and #R = 
1092 for this region. When this value is compared with 
the Reynolds numbers listed in Table V, a glance at 
Figure II shows that the above assumption is valid for 
the waves studied in this thesis. However, if a local 
iR at the bottom of a structure leg becomes less than 
about 100 in ocean waves, the assumption will lead to 
inaccuracies. 

The drag and inertia components were of the same 
order of magnitude for the experimental runs of this 
thesis. Harleman and Shapiro | 27] have found that for 
this particular situation agreement is not good between 
the proposed theoretical forces and the experimental 
forces. Because the theoretical and experimental inertia 
forces of this thesis agree quite closely, (Tables VI 
and IX), it is suspected that the discrepency is due to 
the method of evaluating the drag force component. Ex- 
periments conducted by Wiegel, Beebe, and Moon [ Sal to 
determine average values of Cy in a wave train show 
much scatter, with Cy varying between 0.2 and 2.5 for 
Reynolds numbers of around 10". They suggest that a 
factor of between 2.0 and 2.5 be applied to the drag 


force if an average value of Cy is used. 
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Keulegan and Carpenter /32/, in their work with 


horizontal cylinders, have related C.. to the period 


D 
parameter Ynax? rather than to some Reynolds number. 
If the drag ea is to be averaged over the 
height of a vertical cylinder, the authors suggest that 
Cy and (Cre sO ties be included with A and B in the com- 
puter program (Appendix D). A look-up table for drag 
coefficient, based upon the curve of Cp vs Unax? 


D 
in [32], would have to be included in the program. 


Free Oscillations and Natural Frequencies 
of the Structure Model 


In Appendix H it will be seen that the free os- 
cillation data is somewhat odd in that when the structure 
swings from a negative peak to a positive peak the ampli- 
tude of displacement often increases. Theoretically for 
the damped system this amplitude must decrease [22]. 

When the authors saw this result they immediately sus- 
pected unequal spring constants in the two directions, 
but a test showed that the spring constants were the 
same. The authors are unable to explain this oddity. 

It is seen in Appendix F (or C) that in all con- 
figurations the natural frequency of the model is higher 


in still water than in air. One expects just the reverse, 
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because of viscous damping by the water and the added 
mass effect. The authors cannot explain this result, 
but they conclude that the added mass effect is 
negligible. 

It is also seen that for all runs with the added 
weights the frequency of the braced structure model is 
higher than that for the unbraced structure model , both 
in air and in water. Theoretically spring constant k 
was not changed when the model was braced. Consequently, 
the addition yee bracing mass should have iowered the 
natural frequency Ww, ( = ¥k/m). Because the frequency 
is higher, however, for the braced model, it is probable 
that the bracing increased the stiffness k slightly. The 
authors did not suspect this result at the time of the 
experiment, however, and consequently stiffness k was 


not measured for the braced model. 


2. Computer 


It is mentioned in Appendix D that velocity profiles 
at crest and trough regions were written as output for 
the computer program used to develop Figures III, IV and 
V. These profiles, most of which are not submitted in 
this thesis (see comment, p. 103), show some interesting 


aspects relative to the boundary limits on Figure V. The 
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profiles of particular interest are those for the 
trough region. The computer output gives values of 


Uu Se 


trough We trough 


Vea : 


For certain critical values of wave parameters 
g/L and H/L, ~~ peak value of Ptrough| is shown to 
occur at some level other than at the surface, some- 
times at the ocean bottom. In the extreme situation 


at the surface (= U ) becomes positive. 


“trough trough 
Table XI gives the values of Ys rough associated with 
these critical parameters. nwo €*e ay nenationless 
velocity UL rough refers to the velocity at the surface, 
and Ved YDottom refers to the velocity 


at the ocean bottom y ed under the wave trough. 


The dimensionless parameter (> trough) refers to 
a 


max Se! 
the surface elevation of the trough, and (—EEOUED ) | mit 


refers to the elevation at which the absolute velocity 
is either a maximum or zero, aS appropriate. Values of 
ad/H and L/H are included in Table XI because the com- 
puter output is identified by these terms. 
The authors suggest that the reader refer to 
Figure V while interpreting the information in Table XI. 
First, for a/L = 0.07 and H/L = 0.04, it is seen 
from Table XI that the predicted height of the trough is 


above the still water level. This result, then, can be 
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B71, 


0.20 
0.25 
0.10 
e15 
2.0/7 
0.08 
e210 


TABLE XI - Trough Velocities Associated with Critical 
Wave Parameters 


O 
O 
0 
0.03 
0 
0 
6) 


L/H 


ORO 
7.14 
Pe 
i250 
2 ooo 
20.00 
22 50 


Verough 


Tan 


1809 


oD 
.0960 
51 5iste' 
.1269 
0788 
.0685 


(_ trough) — (trough) 
Ot -0.1636 680 
0.739 -0.1592 650 

-~ -0.1057 -- 
~- -0.1904 -- 
(630) +0.0021 O50 
0.976 -0.0646 .650 
0.901 -0.1686 2175 


discarded as being outside the range of possibility. It 


is interesting to note that this situation is reached when 


Stokes' 3rd approximation predicts positive velocities in 


the trough region from surface to bottom. 


If the remainder of the data is compared with the 


curve families of Figure V, the following observations 


are noted: 


Each curve of the curve familtes reaches a peak 


value, which occurs approximately where Stokes' 


3rd approximation predicts a positive velocity 


S! 
for trough regions above some (Se) ot: 
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Values of the parameters corresponding to the 
above situation are seen to fall approximately 

on line H/d = 0.653, which has been stated to 

be the theoretical limit for wave characteristics 
in the ocean. The profile of positive velocity, 
beginning at the surface and extending down as 
d/L is decreased, is a result of the increasing 
Significance of mass transport in Stokes! 3rd 


approximation. 


Values of the parameters for which u Ls 
trough 
never positive (and for which its peak value 
does not occur at the surface) are seen to fall 
approximately on the line for which H/d = 0.603. 
This line corresponds to the limit of applica- 


bility of Stokes' theory. 


With regard to the statements made above, it is 


to be noted that, when plotted on the curve families 


of Figures III, IV, and V, lines of constant H/d converge 


for low values of d/L. Also, in explanation of ob- 


servation number 1 above, it must be remembered that in 


solving for the force multiplier B from (15) it is 
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a 
S ! 


E ; trougn 
assumed that Ue rough is negative for all — , 


If the proposed procedure is to be extended to regions 


fom which this is not true, then 


2 
trougn 


Yt rough| “trough 


should be substituted for u gle LS) Ie 
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H. Original Experimental Data 


Water temperatures (°F) for the various runs 


are given in Table XII. 


TABLE XII. Water Temperatures (°F) 


Wave en = Sabie Cyl inate 
A 62 Sie 63 
B 62 61 63 
C 62 61 63 
D 62 61 op 


A typical wave gage calibration curve is shown in 
Figure XVIII. A typical force gage calibration curve 
is shown in Figure XIX. "Record reading" on these curves 
refers to the deflection of the stylus needle on the San- 
born recorder. The wave gage is calibrated about the 
Still water level, while the force gage is calibrated 
about the zero force position. Typical Sanborn recorder 
wave profile traces are shown in Figures XX and XXI. 

The damped natural frequency of the 0.5 in. dia. 
cylinder/force-gage combination was 4.8 cps, while for 
the 1 in. dia. cylinder/force-gage combination (used 


for transverse forces) it was 17.85 cps, 


152 





Figure XVIII 
. Typical Wave Gage Calibration 
(for B Wave Profile) 
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Properties of the Structure Model 


Experimental data for free (natural) oscillations 
of the model is given in Table ALI1L. Run 3, for example, 
designates a run in air with the unbraced model with 3 
pounds of added weight on the platform. The letter Z 
designates the braced model. The letter W indicates runs 
with the model in still water. The tabulated numbers are 
peak displacements as obtained from movies. A number may 
not represent the actual peak, which might have occurred 
between frame exposures. After viewing the movies the 
authors deduced that in those few cases in which an error 
could exist it could be only about 0.01 inch. The elapsed 
time listed is between the first and last recorded peak- 
displacements. 

The weight of plastic in the unbraced model (in- 
cluding base) is 22-3/4 lb. The weight of plastic bracing 
is about 0.4 1b. The weight of the bolt is 0.200 lb. 

The total length of bracing is 186 inches. The 
length of bracing Y, below the still water level is 
149 inches. 
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TABLE XIII - 


Run 
Time: pyle! 
Peak -.89 
Displ 98 
-,68 
»(O 
- 48 
RSul 
-.36 
Run OZ 
Time 3.16 
Peak -.67 
Displ. 74 
- 48 
54 
-.36 
41 
-.25 


Experimental Data for Free (Natural ) 


Oscillations of the Structure Uscitiations of the Structure Model 
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Il. Vortex Shedding on a Vertical Cylinder in a Wave Train 


General 


It has been mentioned in this thesis that 
in the development of a practical offshore structure de- 
Sign transverse or "lift" forces caused by vortex shedding 
may be included. It is to be expected that these reactive 
forces are periodic in nature, with magnitudes and fre- 
quencies being functionally related to the horizontal com- 
ponent of wave orbital particle velocity u. Furthermore, 
Since u = 0 at discrete values of wave phase angle 6, it 
may also be expected that this periodicity in vortex 
shedding is repetitive with the passage of each wave in 
a regular train. 

As an illustration let 
F,. = magnitude of transverse force caused by a vortex 
being shed 


O55 = wave phase angle corresponding to Foy 


where i=1,2,...nmn = sequential number of a wave 
in a regular wave train (i = 1 
for the first wave, etc.) 


j= 1,2,...m = sequential number of the vortex 
shed during the passage of one 
wave (j = 1 for the first vortex 
shed during the passage of each 
wave, etc. 
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Then 
F,. = constant (i = 1,n) 


9 sce constant (i = 1,n) 


but since P55 is a function of Uy 5? 


Fi, 7% constant (j = 1,m) 


Thus it is seen that, because of the repetitive 
nature of these forces, if their magnitudes and fre- 
quencies can be predicted a Fourier series representation 
could be utilized in developing a procedure for calcu- 
lating transverse platform displacements. This procedure 
would not differ appreciably from that presented in this 
thesis for longitudinal wave forces. It is the intent of 
this appendix to present an approach for the determination 
of these vortex forces and snedding frequencies and to 
show the results of experiments performed by the authors 


to substantiate this approach. 


A, Theory 


Transverse Forces 


Before the transverse forces are discussed, 
it might be well first to make some qualitative remarks 
concerning the significance of considering these forces 


in a structure design. Mention has been made (Introduction) 
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of the reported failure, due to fatigue, of a piling sub- 
jected to transverse vibrations in ocean waves. 

It is to be noted also that the frequency of this 
force will be at least equal to wave frequency, and for 
very long ocean waves it could be many times greater 
than wave frequency. It is possible that this force 
could create a resonance situation in an offshore 
structure. 

Nothing has been found in the literature relative 
to the magnitudes of tnese forces as a vertical cylinder 
in a wave train, and, to the authors! knowledge, no 
measurements of the forces have been made in surface 
waves. However, if one relates the magnitude of this 
force to the maximum drag force by applying other type 
measurements and steady state theory, some interesting 


observations can be made. 


Let 
F, = maximum transverse force acting 
on vertical cylinder of length 
Fy = maximum longitudinal drag force acting 
on vertical cylinder of length 
Then 3 
U 
dF, = C f D “max af (66) 
2 
=c.Zip2 
dF, = C, @D max dd (67) 
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In steady state flow Landweber [ 33] deduced that 

Cr, = 3.44 3 4 Cy» while for a vertical cylinder moving 
longitudinally with Simple harmonic motion in a stationary 
liquid Laird | 34] found Cy = 0.357 by measurement. lLaird's 
measurement was for a cylinder diameter of 1" and with 
Umax ~ 17 ft./sec. Laird, Johnson, and Walker / 35] 
found experimentally that C, - 0.3 ¢) for a 2" diameter 
cylinder moving under the same circumstances as mentioned 


above. If C, is taken as approximately equal to 1, this 


D 
observation substantiates the value of Cy =z 0.37. The 


above values of C, represent the extremes which the 


L 
literature gives as bases for comparing transverse and 
drag forces. 

It is suspected that the value of C. given by 
Laird most closely approximates the real situation. 
However, two important factors must be considered. First, 
the artifice employed by Laird in approximating a vertical 
cylinder in surface waves is deficient in that the water 
velocity past his cylinder is not a function of the distance 
from the water surface. In the real situation not only is 
Umax €ontinually varying from bottom to oe but also, 
as a wave crest passes a cylinder, the length over nition 


any force acts is not the water depth d but rather d +o 


(%o = height of wave crest above the still water level). 


141 





Secondly, the experimental value of Cy = 0.37 was for a 
given cylinder diameter D and given eon” It is to be 


expected that C like C is a variable with dependence 


L’ De 
on some defined Reynolds number. 

From (66) it is seen that dF, = p(aes From the 
definition of Us it is seen then that dF, is that in- 
cremental force which would result if a vortex were shed 
at the instant that the crest passed the cylinder. Gen- 
erally the vortex is not shed exactly under the crest. 
However, if one is interested in the maximum force which 
might occur, as is the designer, he can justify an as- 
Sumption that the vortex is shed under the crest, where 
the particle horizontal velocities are maximum. There- 
fore, it is assumed that a vortex is shed when 9 = 90°. 

Because waves of finite height are being considered, 
it must be noted further that the design force for the 
crest region will not be equal in magnitude to the design 
force for the trough region. 

As a further complication it was noted by the authors 
in the experimental phase of this thesis that some 
vortices shed in the crest region were swept back over the 
cylinder by velocity reversal during trough passage. These 


cause significant additional forces on the cylinder, but 


how they add, detract, or interfere with new vortices about 
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to be shed is a matter of conjecture. By the time these 
vortices arrive back at the cylinder, some of their vorticity 
has been diffused out of the vortex street. Thus it is 
difficult to determine if the reduction in force in the 
trougn region is the result of new vortices being shed 
in a region of lower orbital velocity or if the re- 
duction is caused by vorticity diffusion of the swept- 
back vortices. Visual observation must be made at,the 
same instant that force is recorded to resolve this issue. 

In the ocean, where for long waves u is in the 
same direction for a relatively long period of time, it 
is proposed that most forces resulting in the wave trough 
will be caused by new vortices being shed, with only re- 
latively few crest-shed vortices being swept back by 
velocity reversal - these few being immediately after 
9 = 180°. The circumstances imposed by laboratory 
techniques where short wave periods are used make it im- 
practical to verify predications on transverse forces 
in way of wave troughs. The authors feel that field tests 
mst be employed for .this type of verification. 

For this reason, it is suggested that if a value 


for F, for the wave crest region is predicted and 


L 
verified experimentally, then a value for tne trough region 


equal in magnitude to K F, (K &—1) could be used in the 
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design stage. It is proposed further that for the 
Structure displacement solution a forcing function 


of the following type be used. 


i= © < 180° 


F= F, sinw, t (68 ) 


tee = 9 4 360° 


Fe KF, sinw, t (69) 


where og vortex shedding circular frequency in 
crest region. 


wo, = vortex shedding circular frequency in 
trough region. | 


K = constant (K < 1) whose value is dependent 
upon the assumed origin of the force. 


More will be said about shedding frequencies later. 
The sine function results from the assumption [22] that 
a vortex shedding force i8 varying sinusoidally. 

The remainder of this section on forces will be 


concerned with the prediction of F, and K, where K will 


L 
be that value based upon the assumption of new vortices 
being shed in the trough region. Several references may 


be cited (for example [ 35,94, 50201 | which can be used to 
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pny 
derive values of Cy by the steady state Von Karman 
vortex street analysis, and this procedure will not be 
repeated here. Rather, the previously specified values 


of C, = 3.44 and C. = 0.37 will be used to place limits 


L 
on the predicted force in formula (66). C, = 3.44 must 
be considered to be very conservative, however, since 
its value is based upon the assumption that u = Unax 
for 0° € 9 €< 180°. 


From (66), F, ={C Spu 2 ae (70) 


where the limits of integration are as yet un- 
specified. Before attempting to solve this equation 


for F it might be well to first consider what criterion 


’ 
may be used to indicate the boundaries for which vortex 
snedding forces are to be considered. For steady flow 
a regular Karman vortex street in the wake of cylinders 
has been observed only in the range of Reynolds numbers 
IR (= =) from about 60 to 5000. At lower values of R 
the flow is laminar, and no vortices are shed. For 
higher IR there is complete turbulent mixing [18]. 

For the non-steady flow associated with wave 
passage, Keulegan and Carpenter (32 / have developed a 
more realistic dimensionless parameter called the period 


parameter. From experiments with horizontal cylinders 


located at various depths in a sinusoidal wave train, 
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they found that the lower boundary for vortex street 
development was coincident with a value of the 


"period parameter" between 12.5 and 15. 


Uu £ 
where period parameter = aE (71) 
u — horizontal particle velocity under 


aa crest at the cylinder depth 


T = wave period 


The significance of the period parameter lies in 
Pirerttical distance a particle must travel in one 
direction for a given velocity in order to produce a 
vortex. For a complete discussion of this approach 


the reader is referred to the original work [32]. 
u 





In this thesis the conservative value of = = 12.5 
will be used for the lower boundary of vortex shedding, 
while it is noted that Keulegen and Carpenter observed 
regular shedding for na as high as 110. 

Two alternatives are therefore proposed for supplying 
the limits of integration to (70). First, it may be con- 
sidered that the vortex extends from the surface down only 
to the depth where Snax” = 12.5. Secondly, it may be con- 
sidered that the vortex, once started on the surface, 
extends all the way to the bottom of me cylinder. This 


u 
second proposal presupposes that = = 12.5 is the 





critical value for starting a vortex on the surface, but, 


146 


once the starting inertia is overcome, the period para- 
meter is no longer critical in cutting off the vortex 


at a specified point along the cylinder. 
T 


U 
Regardless of which proposal is followed, as 2D 


aga Ue 2. at surface) affords the limiting boundary 
below which no vortices are shed. Intuitively the 
authors feel that the true solution is a compromise be- 
tween the two proposals outlined above. However, both 
proposals will be tried for the specific waves of this 
thesis. 

For descriptive purposes the solutions proposed 


thus far are designated as follows: 


(a) vortex extends to cylinder bottom. 


a +%o 


5 Pg) 
F, = Ci 5 Du, ad for crest region (72) 
y 
d-(H- %o) 
ae 
e Ca eal af for trough 
L L 2 trough region (73) 


where y is the distance of the bottom of the cylinder 


above the bottom of the tank (or ocean). 
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AL 
(b) vortex extends only to a depth where — = 1235 


= 12.5 for crest and trough regions, respectively. 





ere mil 
bad pough 


d+%o 
Jae 
= J Cy 5 D max af for crest region (74) 
vo 
a-(H-%o) f 
i eee off 2 b 
P, f C15 D UEP olan a4’ for trough region (75) 
sf 


Uu eg 
where y' is that depth where — = 12.5 and 


u sill 
—ErOUEn 12.5 for (74) and (75), respectively. 


(1) ¢ 


_— 3.44 (Landweber - steady flow) (76) 


(2) C, = 0.37 (Laird - simple harmonic motion) (77) 


As an example, 
solution (a-2) refers to (72) and (73), using a value 
of C, from (77). As mentioned previously, it is expected 
that the true solution lies between (a-2) and (b-2) with 
solutions (a-1) and (b-1) being the conservative extremes. 
The solutions of (72,73,74,75) are of the following 


form: 
PL = cop SP (78) 
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waere 

The is tne mean-square velocity over tne Lengen to 
be considered. If the vortex and the cylinder extend 
to the bottom of the water u° can be determined directly 


from (20) and (21). 
(20) 
(21) 


where A is determined from Figure IV and B, from Figure V. 
inh (2d 
fo, 2 (1 +3 fq Se GH 
at 5 O/H (sinh J) (6) 


The parameter “C2 nay also be read from the computer out- 
put (not plotted), Appendix D. If the vortex does not 
extend to the bottom the following simplification may be 
Meea, heirerence [17] Shows that for distances removed 
from the bottom ue vs. S is approximately linear, where S 
is the distance above the bottom. A conservative value 
of uc may then be taken as 


oe 
a? . Ui tu (utd = 12.5) (79) 


a 
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This equation requires use of (10) and (11) or the com- 
puter output in Appendix D only two times to get a. 
The authors feel that any loss of accuracy is justified 
when tne procedure as a whole is viewed. 

The proposed solutions therefore reduce to the 


following forms: 


(a-1) crest: F, = 3.44 Z D Wax] 4 + Le - a (80) 
trough: F, = 3.44 fp | d - H -'"7,) - y] (eae) 
(a-2) crest: F, = 0.37 Zp "a [a + <i - y ] (82) 
trough: F, = 0.37 SD Veron fa - (H-%) - 9] (83) 
? crest: F, = 3.44 Z D Trax | + We = yt] (84 ) 
trough: F, = 3,4 a D Ue rough) a (H - Te? - y'] (85) 
(b-2) crest: F, = 0.37 £ D a x) 4 = y' | (36 ) 
trough: F, = 0.37% D Te rough] 2 - (H =) ~ y'] (87) 
y = distance of bottom of cylinder above water 
bottom 


y' = distance from bottom to point where = =x 12.5 
and the u“'s refer to the mean square velocity over the 


vortex length. 
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The. force centroids are exactly the same as the drag 
force centroids (Figure III) if the vortex extends to the 


bottom,and for any other vortex length the centroids may 
2 


be assumed in a manner analogous to the assumed Wa and 
er oucn” 
Shedding Frequencies 
Let 5. = cylinder Strouhal number for steady flow 
SS. = Ys D (88 ) 


Cc cw ou 


where W is the vortex shedding circular frequency and 

u is the steady state velocihty. For values of IR between 
10? and 10" experiments by Roshko [26] Show ancaverage 
value of S, of about 0.206. This is the range of concern 
in the experimental work of this thesis. For other ranges 
of Reynolds numbers, one is referred to the reference 
cited above. 


Using 5S, = 0.206, one obtains for tu.» 


0.206 u (29) 
tk (89) 
From (89) it is seen that in a wave train the 
shedding frequency (and number of vortices shed per wave 


cycle) is a function of the orbital particle velocity on 


-~ 
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the: surface of the water (U). This statement is made 

because vortex shedding is initiated at the surface. 

The question next arises as to what value to assign to U. 
Laird [34] found that for his cylinders moving 

in simple harmonic motion no vortices were shed until 

the velocity reached 70% of the peak velocity. If this 

result is to be considered it seems reasonable that (89) 


may be written as 


U Lege U 
Ww, de - + “max) (er) _ EE max) (27) (99) 


(0.7 Vax “ Uinax ) 
where : 5 = "average velocity" at the 


surface during vortex 
shedding in the crest 
region 


and 


0.206(0.85 “Ytrougn) (an) 
Ww, ™ D = 


(91) 
for the trough region. | 


Keulegan and Carpenter i pe); on the other hand, in their 
work on horizontal cylinders in sinusoidal waves propose 
that the shedding frequency be related to the average 
velocity over a half cycle. If tnis reasoning is applied 
to (89) the result is 


U 
Wo 3 pa) a (92) 


for the crest 
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and 


U 
W, i Le? (93) 


for the trough 


Tne authors feel that it would be reasonable to use 

(90) and (91) for laboratory work and for steep waves 
in the ocean, while (92) and (93) may be used for long, 
low ocean waves where the period is long enough to per- 
mit many vortices to be shed. Equations (90) and (91) 
are used by the authors to obtain shedding frequency. 
The designer may want to modify this saqrete Recee however, 
when it is realized that the worst frequency to use is 
that frequency which most closely approaches the natural 
frequency W,, of the structure. In any case the vortex 
force, of this frequency, would be applied to the legs 
of tne structure only for the time duration between a 
being equal to 12.5 on either side of the crest and 
trough. The forcing function proposed by the authors 


assumes the force to be active for the entire wave 


passage. 
Let (i). = time required to shed 1 vortex in 
the crest = 27 (94) 
eu, , 


(Ts), = time required to shed 1 vortex in the 
: trough = 21 (95) 
2), 
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Then from (90) and (91) 


: 86 D 
; Sc 2(0.206)(0.85)U_. ver | 
ond 

on a (97) 


s°t Uerough 


= = time for passage of wave crest or wave trough. 


n= number of vortices shed during crest passage. 


number of vortices shed during trough passage. 


U 


= r = z Max 
me “SRE, S 7D (98) 


U 
2d t “trough 
net ler, = SED (99) 


Because n. and n, must be integral numbers, it is pro- 
posed to accept the next lower whole numbers in lieu 
of the values computed from (98) and (99). These inte- 


gral values of n, and n, then give 


Dy, ag (100) 


where W is th2 wave circular frequency. 
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The proposed theoretical procedure is to use (80) 
through (87) to determine transverse forces, and (98) 


through (101) to determine the forcing frequencies. 


B. Results, Theoretical and Experimental 


For the experimental runs the same waves (A, B, 
C, and D) were used that were used for longitudinal 
force determination. Transverse force measurements 
were attempted by using the 0.5 in. aie cylinder 
described in Section II-B. However, it was found that 
the B, C, and D waves caused resonant vibrations of the 
cylinder/force-gage combination in the transverse 
direction, with the A wave being the only wave which 
gave usable results. For comparative purposes the 
transverse forces caused by the A wave are plotted with 
the longitudinal forces and the results are shown in 
Figure xVila. 

As a consequence of these results it was decided 
to center investigation on alin. dia. cylinder. 
This cylinder cleared the bottom of the tahk by 
12-5/8". Therefore, in (80,81,82,83) y is taken as 
12-5/8" = 1.052 ft. By measurement the damped natural 


frequency of the 1 in. dia. cylinder arrangement was 
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17.85 cps, as compared to 4.8 cps for the 0.5 in. dia. 
cylinder arrangement. 

Before the results are presented it would be in- 
teresting to predict whether or not vortex shedding is 


expected for the various runs. Table XIV lists values 





Cee U fc 
of wave period, as well as = and rpuen for 
the 4 waves of the experiment. Wave periods are de- 
termined from (8), while the velocities are taken 
directly from the computer output (Appendix D). All 
other values in the following analysis are taken from 


Appendix D, unless a particular equation is specified. 


TABLE xiV- Wave Periods and Period Parameters 











Wave T (sec. ) “nax ~trougi 
A eeULs 15 ee er ak rae Eas, 
B 1.558 15.45 -12.00 
C Lie. DSO dihaih -~12.53 
D 1.200 18.90 -12.90 


Since the critical value of = is considered 


to be 12.5, it is doubtful that the trough region 


can support vortex shedding. 


In faet, all forces in 


the trough region were observed to be caused (in part 


at least) by crest-shed vortices being swept back to 
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the cylinder during velocity reversal. 
shed vortices are not considered further, except to 
State that for the runs performed by the authors, 
we WD, - | 

Parameters used in the solutinn of (80, 82, 84, 
86) are shown in tabulated form in Table XV. The 


parameter ( refers to the velocity for which 


ree’ bet 
the period parameter is equal to 12.5 and is used as 


the entering argument in Appendix D to determine y'. 


The parameter ey is the velocity under the crest 


at the bottom of the cylinder. The parameter (a,.°), 


is for use with (80) and (82) and is equal to 
2 4 


peak + (“max by = 
D .  #The parameter (Clemo is for use 
with (84) and (86) and is equal to 
Z 2 
y « crit 
oS —— . Velocities are in ft./sec. 


Hence, trough- 


TABLE XV ~- Wave Parameters for Transverse Force Determination 


mae Mo (Ft) Onax (Unaxorit Ca ax! et) Tne 
A 0.12 0.568 0.518 0.399 2.09 0.241 
B ery 0,025 0.669 0.452 2,00 0.4402 
C 0.200 1.050 0.774 O62. 1.97 § 01657 
D Omee7 12510 0.870 0.453 1.95 0.958 
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The transverse forces (F, ) as determined experi- 
mentally are shown plotted with the corresponding wave 
profiles in Figures XXII a-d. Table XVI compares 
theoretical and experimental results for values of FL 
and n, (number in parenthesis is integral value). The 
value of w, (= uw, for reasons mentioned) is deter- 
mined from (100). The theoretical values of n. are 
determined from (98), and the various F,'s from the 
appropriate formulas (80, 82, 84, 86), with inputs 
from Table V. 


f = 1.94 slugs/ft.? d = 2.25 ft. 


TABLE XVI - Comparison of Theoretical and Experimental 
Transverse Force Results 





Fy un/em (cps ) 


Hon. | 
U & ath roe ar | 
Method Theory Exp .| Theory Exp. | 
pare l| a-1 | a-2 | b-1 b-2 |Exp. 


A 110.0883 !0.0095 |0.02350;0.00248] 0.014)//2.4(2) 2 0.993 10.993 








B |j0.1680 jO.0181 |0. 0651110, 00704 0.072H2.7(2) 2]/ 1.283 |1.283 
C |]0.2550 jo. 027 iO. 1167/0. 01255] 0.12813.0(2] 2]| 2.220 11.482 |, 
D 


0.3800 jO.0408 '0.1950 [0.02100] 0.088}|3.3(2] 2|| 2.500 f1.667 


C. Discussion of Results 


From Table XVI and Figures XXII a-d it is seen 
that the assumption of vortex shedding frequency 
being equal to an integral multiple of wave frequency is 


Wye 
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valid. Laird [34 } found experimentally that vortex 
shedding appears to occur in pairs, and this may be 
why the experimental values of n, for the C and D waves 
are equal to 2 rather than 35. For all experimental 
runs trough shedding frequency wW, was equal to crest 
shedding frequency W 2 as explained previously. 

The vortex forces in the trougn region were not 
determined analytically, because, for the experimental 
runs of this thesis, these forces are caused (in part) 
by back-sweeping crest-shed vortices from which some 
of the vorticity has been diffused out of the street. 
From Figure XVilaforthe 0.5 in. dia. cylinder, it is 
seen that the experimental maximum vortex force is 
43.4% of the experimental longitudinal total force and 
90% of the experimental longitudinal drag force. This 
90% value is considerably higher than the 30% value 
found experimentally by Laird, Johnson, and Walker 25. 

From Table XVI it is seen that method (b-1) agrees 
most closely with the experimental force results. Method 
(b-1) refers to (84), where C; is taken as 3.44, and 
the vortex length terminates where Unaxt = 12.5. Methods 

— 
(a-2) and (b-2), which assume C, = 0.37, are seen to give 
forces which are far too small. It is seen, therefore, 
that a method of solution between (a) and (b) together 
with a method between (1) and (2), would give analytical 
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results most closely approximating the experimental re- 
sults of this thesis. The 90% ratio of C,/Cy; determined 
experimentally for the 0.5 in. dia. cylinder (Figure XVIIa), 


~ 


may indicate a basis for selection of C More experi- 


L’ 
mental tests are required to narrow the boundaries for 
the analytical forces, as shown in Table xVI, but 
it is reassuring that the experimental forces lie 


within the range defined by the proposed theories. 


D. Conclusions and Recommendations 


1. The proposed method of determining a design vortex 
shedding frequency is valid. It is suggested that 
the designer pick a value of n, Which is on the 
Side of the computed value to cause WwW, (or Ww) 


to approach the structure natural frequency BY) 


2. Until a more accurate prediction may be made of 
the vortex shedding forces, a safe design is to 


assume Cy. = Cp: 


D. It is recommended that field tests be conducted 
with ocean waves of long period to verify the 
authors' procedure for the case of trough-shed 


vortices. 
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4, If platform rotation is considered, vortex 


shedding forces will assume added importance. 
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